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EXPERIMENTAL INVESTIGATIONS OF OPTICAL 
PROPAGATION IN ATMOSPHERIC TURBULENCE* 

by 

Michael Williams Fitzmaurice 
Goddard Space Flight Center 

INTRODUCTION 

The study of optical propagation in random media has been renewed with 

1 2 

increased vigor in the past two decades. The works of Chernov and Tatarski 
in the Soviet Union during the 1950's rekindled interest and advanced the state of 
knowledge so that many previously unsolved problems became rigorously 
amenable. It appears that these efforts were motivated primarily by a desire 
to fill theoretical voids in the classical treatment of certain problems in electro- 
magnetic theory and acoustics. Following the initial realizations of the laser, 
however, these studies gained an increased utility. It became immediately 
apparent that attainment of the exciting goals forecast for the laser would require 
that it retain it's spatial coherence, temporal coherence, and high energy density 
over the propagation path. Clearly, transmission through the earth's atmosphere 
would introduce a randomness to each of these factors and thereby impose certain 
limits on performance capabilities. During the same time-period, increased 
emphasis on aerial photography from aircraft as well as satellite platforms 

Dissertation submitted to the faculty of the graduate school of the University of Maryland in 
partial fulfillment of the requirements for the degree of Doctor of Philosophy, 1970. 
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generated a need for an understanding of resolution limitations imposed by a 
turbulent atmosphere. In addition, the field of optical astronomy has traditionally 
been hampered by the earth's atmosphere; thus, considerable effort has been ex- 
pended to understand the causes and effects of the turbulent atmosphere upon 

3 “6 

earth-based measurements. For all these reasons, the monographs of 
Chernov and Tatarski aroused considerable interest and have served as a starting 
point for much theoretical and experimental work on the optical propagation 
problem in the decade of the 1960's. In this dissertation we shall — 

(1) Review the causal mechanism for the optical propagation problem (i.e. , 
the random refractive index field). 

(2) Summarize the more important theoretical and experimental results now 
available . 

(3) Present the results of a group of new experiments which, it is hoped, 
will clarify certain aspects of the propagation problem. 

The material is organized such that Chapter 1 deals with items 1 and 2, and 
Chapter 2 presents the new experimental work, item 3. 
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CHAPTER I 


THEORETICAL CONSIDERATIONS 


l.l Optical Propagation Through Turbulence: Problem Definition 

In general, a propagating electromagnetic wave can be perturbed by many 
different mechanisms. Different types of perturbation are treated analytically 
by very distinct techniques which are sometimes exact but often approximate. 

We now review four types of propagation problems so as to define precisely the 
problem with which we are concerned in this paper. 

The first type is the free space propagation problem. Here a wave is launched 
into an evacuated medium and progresses in accordance with the classical wave 
results obtainable from Maxwell's equations. In regions near obstacles, diffrac- 
tion theory with suitable boundary conditions is employed. Usually, Kirchhoff 
diffraction theory^’® is employed, although it is well known that it is an approx- 
imate technique and, in fact, has some internal contradictions. Other 
approaches such as the Rayleigh-Summerfeld theory 9 ’^*' 1 ' 2 are available and 
are sometimes employed. 

In the second type of problem we consider a material medium which is 
characterized by constant properties. That is, the medium has a constant com- 
position, temperature, pressure, and permittivity. No particles with sizes greater 
than the molecular scale are present. In this case, one finds that the molecular 
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matter is excited into oscillation by the incident wave and behaves as a collection 


of classical Hertzian dipoles with radiated fields which are calculable from well- 
known formulas. This type of scattering was initially treated by Lord 
Rayleigh^ 5 >1- (5 an d is widely known as Rayleigh scattering. 

In the third type of problem, we consider a medium where the dominant 
perturbations to a propagating wave are caused by particulate matter with scale 
sizes which are equal to or greater than the wavelength of the incident radiation. 

A problem of this type has been solved by Mie^ for the case of a plane wave 
incident on a spherical particle of arbitrary size and permittivity. The scattered 
energy is expressed in terms of an infinite series of terms representing products 
of associated Legendre polynomials and spherical Bessel functions. Exact 
series and useful and approximate solutions are available.!® In practice the 
usefulness of either result is limited by the knowledge available on particle 
size distribution, density, and composition throughout the propagation medium. 

A great deal of information is available in these areas for particular cases. 20-23 

A fourth type of propagation problem is the subject of this thesis. Here 
the medium is devoid of particulate matter and is characterized by zero con- 
ductivity and certain other average properties (composition, temperature, pres- 
sure, permittivity, etc.) which may be slowly varying functions of position. We 
also allow the point -by -point values of the medium properties to exhibit small 
deviations from the average. 

Physically, this propagation problem may be interpreted as follows : when a 
wave (e. g., plane wave) is launched into the turbulent medium, the small deviations 
in medium properties which result in correspondingly small refractive index 
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fluctuations cause the wavefront to quickly become non-planar. For relatively 
short distances into the medium, the amplitude along constant phase surfaces 
remains essentially constant, although the wavefront distortion grows progres- 
sively larger. It should be clear that, for long enough propagation paths, the 
warping of the equi-phase surfaces will result in a shaping of the amplitude dis- 
tribution across the wave. So, in general, one must solve for the complex-valued 
field quantity u which satisfies the linear stochastic differential equation* 

V 2 u + k 2 n 2 u = 0, (1) 

where the wave number k - 2 n/\ - u/c and n is the refractive index. 

The field u depends on position and time , i.e., n - n(F, t)- However, we 
account for the time dependence of n by letting each n( r , t . ) , for i = 1,2,..., be 

considered a particular realization of the stochastic process {n( r)} . This as- 

< • 

sumption is clearly valid, since the time constant associated with changing spatial 

distributions of n f^lCT 3 seconds) is much greater than typical wave transit times 

24-26 

through the turbulence and the wave period. 

In principle, one could solve the problem (Eq. 1) for each realization n. (r), 
catalogue the results, and by monitoring the behavior of n( r) over the propagation 
path extract the precise time behavior of the field u. This, of course, is not 
practical, since the number of realizations is infinite [n(T. , t)is continuous] 
and the instrumentation required for obtaining instantaneous 
values of n over the entire path is never available. Instead, we approach the 


The justification for using this equation as the starting point is given in Section 1.3- 
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problem from a statistical standpoint and attempt to evaluate certain statistics 
of u (mean, variance, etc.) in terms of the statistics of the refractive index field 
n(7) over the propagation path. Clearly then, we must inspect the index of re- 
fraction fluctuations which occur in the atmosphere, their parameter dependence, 
and their statistical form. This is carried out in the next section. 


1.2 Refractive Index Statistics in the Turbulent Atmosphere 

The index of refraction of air is expressible in terms of the properties of 

27 

the air through the relations 


77. 6P / 0.0075\ , 

1.0 + — j I 1 t — — J 10 , at optical frequencies, (2) 


and 


77 . 6 P , 3.8 „ , 

1.0 + — j x 10 6 + ( 1 - 0. 0195 z) , at radio f requencies, (3) 


where P is atmospheric pressure in millibars, T is absolute temperature in Kelvin 

degrees, A. is wavelength in microns, z is altitude above sea level in thousands of 

28 

feet, and a standard atmosphere is assumed. We restrict our attention to 
the optical case, Eq. 2. In turbulent flow, the pressure fluctuations have been 
shown to be negligible* and it can therefore immediately be concluded that devia- 
tions in the optical index are caused by microthermal fluctuations. Our goal, 
then, is to determine the statistical nature of the thermal fluctuations. 


This point is discussed in detail in Ref. 24, pp. 208-210. 
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The approach used by Tatarski (which he states was initially formulated by 
29 30 

Obukhov and Kolmogorov ) first investigates the velocity fluctuations in a 
highly developed turbulent flow. Highly developed turbulent flow means simply 
that the Reynolds number which characterizes the flow is much greater than the 
Reynolds number at which the transition from laminar to turbulent flow takes 
place. The results are most conveniently expressed in terms of the structure 
function D ik (7) of the velocity field defined as 

D ik (7) - (V.- V.')(V k - V k '), for i,k=l, 2, 3, (4) 

where the V i are the components of the velocity vector in the x, y, z directions 
at the point 7 x , and the V.' are the same quantities at the point 7 x 4 7 . So, in 
general, D ; k (7) is a set of nine functions which describes the intensity of the 
velocity fluctuations of spatial extent r or less. Assuming local isotropy, 
express D. k (7) as a linear combination of D rr ( r) andD tt ( r), where 

D rr (0 = ( v r - v ;) 2 

and 

D tt( r ) = ( v t- v t ') 2 , 

and where V r is the component of the velocity field at 7j in the direction 7, and 
V r ' is the same quantity at the other measurement point 7 t + 7 . Similarly V t is 
the component of the velocity vector at 7j in a direction perpendicular to 7, and 
V t ' is the corresponding quantity at the point 7 t +7. Assuming incompressible 
flow, Tatarski shows through dimensional analysis that both D r r ( r) and D t t ( r) 
grow as r 2/3 for spacings r much greater than the inner scale of turbulence and 
much less than the outer scale. 
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We carry out the analysis of temperature fluctuations by considering 
temperature to be a conservative passive additive to the highly developed turbu- 
lent flow. The problem of determining the statistics of the resulting micro- 
thermal field was solved by Obukhov 3 - 1 - and Yaglom 32 and yielded results for the 
structure function of the temperature field: 


D T (r) = a 2 r 2/3 , for * o «r«L 0 , 

and 

1 N , „ 

D x( r ) = 7 D r ’ for r <<^ 0 ’ 


(5) 


where a is a numerical constant, N is the dissipation rate of temperature fluctua- 
tions due to molecular diffusion, e is the corresponding dissipation rate caused 
by fluid viscosity, and D is the molecular diffusion coefficient. As before, the 
"2/3 law" is seen to hold for separations between the inner [i o ) and outer (L o ) 
scales. The spatial range defined by Z o and L o is called the inertial subrange 
and is considered to denote the limits over which temperature (and, therefore, 
refractive index) fluctuations may be considered isotropic. Proportionality 
constants are usually defined over the inertial subrange as follows: 

D t( 0 = [T(r x ) - T(7 1 +7)] 2 - C T 2 r 2 ' 3 

and (6) 

D n( r > = [ n ( T l) “ n ( T l + 7 )] 2 = C n f2/3 > 

where C T 2 is called the temperature structure constant and C n 2 is called the re- 
fractive index structure constant. The notational change from D(7) to D( r) results 
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from the assumption of isotropy. In general, at visible wavelengths and moderate 
average temperatures, 

C n 2 % 10" 12 C T 2 • 

There are several points worthy of discussion as regards the Obukhov - 
Kolmogorov turbulence model. First, and perhaps most important, there is 
considerable experimental evidence both in the Soviet Union 33 ’ 34 an( j the United 
States, 35 ’ 36 that the model closely agrees with observations in most cases. 
Because of the model itself having a firm basis in statistical turbulence 
theory, it would seem to be very risky to use any significantly different model 
in propagation calculations if agreement between theory and experiment is to be 
expected. 

It must be noted, however, that this model says little about the form of D n ( r) 
for r values near the inner and outer scales. In many problems, the exact be- 
havior of D n ( r) near r = L o does not appear to be very important. 37 ’ 33 How- 
ever, the form of D n ( r) near r = % o may be quite significant 35 and could 
appreciably effect some calculations. Also, the mathematical form of the 

Obukhov - Kolmogorov model appreciably complicates calculations, and solutions 

40 

are obtainable only through computer-processed numerical integrations. 

In an attempt to avoid or alleviate certain of these problems, various the- 
oreticians^’^ ■'■have employed a gaussian model for the covariance function of 
the index of refraction. For example, Beckmann^ 2 used the following covariance 
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model: 


C(r) 


fj. exp 


[- (x ) 2 - at - & 


(7) 


where m 2 is the variance of refractive index field; x, y, and z are the rectangular 
components of r; and X, Y, and Z are the correlation distances in the respective 
directions. Clearly, such a model eases the mathematical complexity in making 
calculations and, in addition, has the decided advantage of permitting anisotropy 
in the refractive index field. However, the model has no basis in turbulence 
theory and has only one scale length per direction. In this sense, the gaussian 
model is far less acceptable than the Kolmogorov model. 

In conclusion, it seems that if one desires to make meaningful comparisons 
between experiment and theory, the observations must be compared with theo- 
retical results calculated through "2/3 law" turbulence models. Any agreement 
of experiment with gaussian -based theory can probably be regarded as fortuitous. 
In the remainder of this paper, we emphasize only those results from propagation 
theory which are based on the Obukhov -Kolmogorov Turbulence Theory. 

1.3 Theoretical Approaches to the Optical Propagation Problem 

Solutions to the general problem of propagation through random media have 
been sought by many theoreticians using very different approaches. We make 
no attempt to review all existent techniques, but rather restrict ourselves to 
those which are most widely used. Each will be reviewed from first principles 


10 



and its limitations noted. Specifically, we concern ourselves with Geometrical 
Optics, Born Theory, and Rytov Theory. 

A common starting point exists for each technique. We consider a non- 
conducting medium without sources which is characterized by a position-dependent 
scalar permittivity, e = e(x, y, z) , and a constant scalar permeability, At = • 

Then Maxwell's equations can be written in MKS units as 


a. V • D = 0, 


c. V • B - 0, 


b. V x E 


dB 
d t ’ 


dD 

d. V x H - -gY ’ 


with the constitutive relations (8) 


e. 

D = 

eE 

f . 

B = 

P-H 


Using the vector identity V x V x A = V(V • A) - V 2 A and taking the curl of Eq. 8b, 
we obtain 


V(V • E) - V 2 E 



( 9 ) 


Assuming harmonic time -dependence (e - ,co 1 ) for all field quantities, and putting 
8d into the right-hand side of Eq. 9, we get 

V(V • E) _ V 2 E = j cojj. ( — jo-'eE) = peE . (10) 
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ilium 111 ii i mi Hi i i 


Using the identity V • (VA) - W • A + A ■ VV applied to Eq. 8a, we get 

V • eE = eV • E + E • Ve . 


Thus, 


V • E 


V • eE - E • Ve 
e 


(ID 


Substituting into Eq, 10, we write 

V 2 E - ai 2 /^dE - ~ • 


( 12 ) 


The index of refraction n for a medium is defined in terms of a phase velocity 
ratio as 


n 


c 

Y 




where u 0 , , and c are the free space values and we assume^ - u 0 . Setting 

the permeability equal to unity and introducing the free space wave number 

OJ 2n 

k = — = — ' 


we can write Eq. 12 as 


V 2 E + k 2 n 2 E = - 2V(E • V In n) 


(13) 
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This is the general equation which must be solved for propagation in a medium 
characterized by a random -point function for the refractive index. Eq. 13 is a 
vector equation and hence is really three scalar equations. Consider the x com- 
ponent and expand the right-hand side. Then, 


V 2 E + k 


3E 


2 n 2 E 


= E. 


dx 


2 /« n + ^ x 


x d In n 


d x 


+ E 


a 2 

dxdy 


In n 


(9E q ^ 2 

~dx dy ln n + E z dxdy ln n 


d In n 
dx d z 


(14) 


Clearly the term V(E • V In n) in Eq. 13 couples the various components of the 
field together and is responsible for polarization fluctuations. Usually this term 
is dropped, and polarization fluctuations are considered negligible on the basis 
of the following argument. 43 Since the deviations of n from its mean value are 
very small (typically, ty^n) ~10 -6 ), then one can assume that the deviation of E 
from the value E 0 it would assume in a turbulence -free region is also small. 

In this case, a perturbation-type solution can be sought for E, i.e. , 

E = E 0 + Ej , 

with (15) 

n - (n> + nj . 

Substituting Eq. 15 into Eq. 13, expanding, simplifying, and dropping all terms 
greater than first order, we find that 

V 2 Ej + k 2 Ej + 2k 2 rij E 0 + 2V(I 0 'V ni ) = 0 . (16) 
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Comparing the third and fourth terms in Eq. 16, we find that the third divided 
by the fourth is approximately equal to ^ 0/ A- Since the assumption -l Q » k ap- 
pears reasonable for all optical problems, the fourth term in Eq. 16 can be 
dropped. Since this term was generated by the right-hand side of Eq. 13, then 
V(E • V In n) in Eq. 13 can be neglected. We conclude that the equation to be 
solved is 

V 2 E + k 2 n 2 E = 0 . ( 17) 

The solution of this equation using the three previously listed techniques is the 
subject of the following sections. 

a. Geometrical Optics . We first note that Eq. 17 is equivalent to three 
identical scalar equations, and so we need only deal with 

V 2 E J ' k 2 n 2 E = 0 , (18) 

where E is any of the rectangular components of E. Eq. 18 can be rewritten in 
7 

the form 

V 2 E + k 2 n 2 E = '—Y L + k 2 n 2 = V 2 In E + | V In E| 2 + k 2 n 2 = 0, (18A) 

where the various forms are easily equated by direct expansion. Let 

E(r) = A(r) exp [i S( r)] . 


Then we have 

V 2 {In A+iS) + |V(/«A+iS)| 2 + k 2 n 2 = 0 
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Separate the above equation into real and imaginary parts and equate both to 
zero. For the real part, 

V 2 /« A + | V In A| 2 ~ | Vs| 2 + k 2 n 2 = 0 . 

For the imaginary part, 

V 2 S + 2 V In A • VS = 0 . (19) 

Noting the relationships used above, we can write 

V 2 A 

| Vs! 2 - k 2 n 2 + —fi— > (20) 

but, since the wave amplitude can be expected to change in space no faster than 
1 /-£ , then 

O 7 


v rv -i- o o o 

A ~ JT << k 2 n 2 - k 2 . 

o 

Therefore, Eq. 20 becomes 

|VS| 2 = k 2 n 2 , (21) 

which is the Eiconal equation. Eqs. 19andEq. 21 are the basis for the geo- 
metrical optics approach. The primary limitation of the approach is that ampli- 
tude changes near boundaries are not properly treated. 44 The neglect of the 
V 2 A/A term is tantamount to ignoring diffraction effects. The amplitude fluctua- 
tions which are predicted by geometrical optics for the turbulent w'ave propagation 
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problem are not in any sense due to diffraction effects but rather are caused by 
the lens effect of inhomogeneities along the path. That is, as one moves along 
the propagation path, the phase field has some curvature associated with it at 
every point. The integral of this curvature along the path results in a focusing 
or de-focusing of the on-axis ray bundle, and hence the on-axis energy density 
either increases or decreases. The predictions for phase fluctuations obtained 
through geometrical optics are in general considered to have a wider range of 
validity than those for amplitude statistics.^ This aspect will be considered in 
more detail in later sections. 

b. Born Theory . As before, we begin with the scalar version of Eq. 17, i.e. , 

V 2 E + k 2 n 2 E = 0 . 

Since the refractive index deviations are assumed small (compared to the mean 
value), a perturbation-type solution of the form 

E = Eq + Ej 

is sought, where E 0 is the field that would exist in the absence of turbulence. 

We let n(T) = 1 + rij (T) and note that V 2 E 0 + k 2 E 0 = 0. Eq. 18 can now be writ- 
ten as 

V 2 E 1 + 2k 2 E 0 + k 2 ni 2 E 0 + k 2 Ej + 2^ k 2 E x + k 2 n 1 2 E 1 ~ 0. 

Since << 1 and E t is of the order n x |E 0 |, we can neglect higher-order terms 
with small error and write 

V 2 Ej + k 2 E x = - 2k 2 nj E Q . 


16 



This equation has the solution 


h CO = 2k 2 nj (r')E 0 (T')g(|t- r ' | ) 

J v 


dv 1 


( 22 ) 


where g{| r - r ' | ) is the free space Green's function defined as 


G(|T- r' |) 


1 exp [ik[ r ~ r' 1] 

477 


r - r 


and the integral is over the scattering volume V. Eq. 22 is the basic result of 
single -scatter Born theory. Its usefulness is seriously compromised by the 
restriction that | E* j/j E 0 1 << 1 .* In fact, it is well known that at visible wave- 
lengths this condition is violated for relatively short ranges through mild 
turbulence. 46 

A multiple scattering approach can be formulated from this theory by a 
series expansions of E in the form 

E = E„ + E t + E 2 + E. , (23) 

where | Ej | ~n t | E._j [ . Substitute Eq. 23 into Eq. 18 with n = 1 + nq , expand, 
and group terms of equal order of smallness. Then we obtain the systems of 


•For phase calculations using Born theory it is required that k 2 ^1 « 1, where -t is a 

characteristic scale for the turbulence and L is the path length. 
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equations 


( V 2 + k 2 ) E 0 = 0, 

(V 2 +k 2 ) Ej = - 2k 2 rij E 0 , 

(24) 


(V 2 + k 2 ) Ej 


= - 2k 2 rij E j_ j - k 2 n 1 2 E j _ 2 , for j > 2 


This set of equations has the solution 


(f ) 


If"' 


+ 2k 2 n, 


Ei- 2 C7' )j 


E i -2 ( r ' ) G (l r - r ' | ) dv', (25) 

for i > 2 . 


We now have a set of integral equations to solve, starting with i = 1. Clearly, 
the range of validity of the technique has been extended, but the mathematical 
complexity is increased so much that the possibility of making meaningful calcu- 
lations is doubtful. The accuracy of the technique versus the number of terms in 
the series does not seem to be well established at this time. 

c. Rytov Theory . This approach has been given wide exposure through the 
work of Chernov and Tatarski and continues to be the subject of much discussion. 
The derivation starts with the transformed version of Eq. 18, i.e. , 


V 2 In E + j V In E| 2 + k 2 n 2 = 0 . (18A) 

Let 

E(r) = A(r) exp [i S(r)] 

and 

n(7) = 1 + nj (7) , 
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and define the complex phase \ p as 


E(r) - exp - exp [/« A(r) + i S(r)] 


Therefore, 


/« E(r) - /« A( r ) + i S( r ) - 0( r) 


(26) 


Substituting this result into Eq. 18A, we can write 


V 2 0 + [ V0j 2 + k 2 (l +nj) 2 = 0 


(27) 


We set 


Therefore, 


0 - + ' 


E - exp [y» 0 + s/'! ] - E 0 exp [ ] > 


(28) 


where E n is the unscattered wave which satisfies 


V 2 E 0 + k 2 E 0 = 0, 


or 


V 2 /»E 0 t |V/«E 0 | 2 * k 2 = 0 
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This last equation can be written as 


v 2 ^ 0 + I V0 O I 2 + k 2 = o . (29) 

Substituting the first of Eqs. 28 into Eq. 27 and subtracting Eq. 29 from the result, 
we obtain 


V 2 ^ + IW'J 2 + 2V^ q • V^ + k 2 (2n 1 + n 1 2 ) - 0, (30) 

which is the nonlinear Ricatti equation. Since rq « 1, k 2 n L 2 is much less than 
2k 2 rij and can be dropped. Similarly, since j Vi p 1 1 « | V p Q j ~ k , then | Vi p | 2 
« 2 Vi p Q • Vi p x and we also neglect the | Vi/^ | 2 term. Therefore, Eq. 30 can be 
written 


V 2 i/ij + 2 Vi/< 0 • Vi/ij + 2k 2 n t - 0 . (31) 

Eq. 31 is the linearized form of Eq. 30 and is the fundamental equation of Rytov 
theory. The dropping of the term ] Wq | 2 in Eq. 30 is called the "Rytov 
Approximation." We now define W in terms of complex phase quantities as 

W = i/ij e^° . (32) 

Substituting for in Eq. 31 using Eq. 32 and after expanding and simplifying, we 
have 


V 2 W - k 2 W = - 2k 2 nj E 0 , 


(33) 
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which has the solution 


W(7) 


1 


2k 2 rij (r')E 0 ( 7 ' ) G (| r - 7' | ) dv\ 


(34) 


and therefore, 

'Pi ( T ) = ~ - - 2k 2 ni (r ')E 0 (r) G(|T -7' |) dv' . (35) 

( r > •'v 

It is noteworthy that exactly the same calculation is needed to obtain W(7) as 
was needed to find the first term in the Born series. It is tempting to conjecture 
that the realms of validity are therefore equal. However, this reasoning can be 
quickly dismissed when one realizes that the nature of the correction term in 
each case is quite different. In Born theory, the correction term is summed with 
the unscattered field to get the total field, while in Rytov theory the correction to 
complex phase (i.e. 7 X ) appears in the total field in a multiplicative fashion. 
Clearly, the range of validity of the two approaches need not be the same. 

The derivation of Rytov theory just completed contains only one noteworthy 
assumption, namely that | j << | Vi/< 0 1 ~ k . This is equivalent to saying that 
the change in the correction to complex phase over distances of a wavelength is 
much less than 2tt radians. If we recall that X << -£ .this condition would seem 

O 1 

very tolerable and we would conclude that Rytov theory has an extremely large 
range of validity. Actually, this does not seem to be the case. Present thinking 
on the relative merits of the three preceding approaches is discussed in the 
following section. 
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1.4 Comparison of Theoretical Approaches: Relative Range of Validity 

Considering geometrical optics first, we noted earlier that this approach 
neglected diffraction effects. In addition, refractive effects on off-axis ray 
bundles will result in caustics and cause a "spill over" of energy into the on-axis 
ray bundle. This mechanism is also neglected. The usual sufficiency condition 
applied for a geometrical optics approach to be valid is that the path length L be 
much less than -A A . For typical near -ground propagation, this results in an L 
value of a few meters. Taylor ’ has questioned the validity of this criterion 
(i.e.,L << £*/\ ) and raises the possibility that the ray optics technique may be 
valid out to several hundred meters. A complicating factor in the discussion is 
that one must define the statistic of interest before assigning a range of validity 
to a technique. 49 It is clear that phase fluctuations are handled in a more adequate 
fashion than amplitude fluctuations by ray optics. This follows since the answer 
one gets for the phase structure function through ray optics is the same as the 
answer obtained through wave optics. 50 

Perhaps the most straightforward experimental method for finding the 
geometrical optic limit for amplitude analysis would be to transmit two widely 
spaced (e.g.,A\ micrometer) optical beams over a variable -length path. Since 
geometrical optics is the high-frequency approximation to the propagation 
phenomena, it has no wave -number dependence in the calculation of amplitude 
fluctuations; the appearance of a statistical difference in the amplitude fluctua- 
tions of the two beams could be interpreted as the break-down point of the ray 
optic amplitude calculation. This type of experiment has not to the author's 
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knowledge been carried out, and we are without any definitive experimental in- 
vestigation of the problem. However, it can be noted on the basis of the results 
presented in Section 2.2 of this report that over a near -ground horizontal path 
under weak turbulence conditions (i.e. , at night) the valid range for amplitude 
calculations is much less than 1.17 KM. The problem is not as perplexing as it 

might seem, since the geometrical optics approach can be shown to be a special 

51 52 

case of the more general wave treatment; ’ therefore, any needed result 
can be obtained using wave optics (assuming the mathematics are tractable). 

Single -scatter Born theory is seriously hampered by the restriction 
| E,/E 0 | << 1. To gain some insight into the magnitude of the problem, we consider 
Protheroe's 33 measurements of stellar irradiance fluctuations as a function of 
zenith angle. Under the weakest scintillation conditions, (i.e. the star is at 
zenith and chromatic averaging is allowed), the normalized irradiance fluctuation 
(standard deviation divided by mean) is about 36%. So clearly, even this rather 
mildly turbulent path is beyond the scope of single -scatter Born theory. Typical 
near -ground paths exceeding a few hundred meters in length induce visible ampli- 
tude fluctuations which are larger than stellar scintillation and are also not 
amenable to this technique. It is possible, however, that for infrared sources 
the magnitude of the wave fluctuations in many cases may be sufficiently small so 
that single-scatter theory is applicable for these ranges. 

Rytov theory contains only one notable assumption; namely, that the correc- 
tion term to the complex phase changes little over distances of a wavelength. As 
noted earlier, this is not a significant restriction for most problems; therefore, 
one is tempted to conclude that the range of validity is great. However, the method 
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has been questioned both in the Soviet Union 5 ^ -55 (as reported by Strohbehn) and 
in the United States . ^ DeWolf 5 ^ argued that the method included geometrical 
optics and single -scatter Born theory and had a range of validity much greater 
than either. Brown 58 viewed the technique as a perturbation solution to the non- 
linear Ricatti equation; after calculating the next higher term in the solution, 
he concluded that the Rytov technique was valid only in situations where the 
single -scatter theory could be applied. Fried 55 interpreted Coulman's data 55 as 
experimental proof that the Rytov approximation gave accurate answers in situa- 
tions which could not be handled by single -scatter or ray optic theory. Various 
other papers have appeared either defending^ ^ or criticizing^ the usefulness 
of the approach. 

Perhaps the most meaningful conclusion that can be drawn in this situation 
is that the problem is not likely to be clearly and decisively resolved through 
analytical efforts. Properly designed experiments should be able to settle the 
issue and in fact, to this writer, it appears that certain aspects of the problem 
are now resolvable. In view of the data presented in this report as well as other 
recent experimental results it seems clear that Rytov theory accurately pre- 
dicts amplitude statistics in cases where neither geometrical optics nor single- 
scatter theory is applicable. The remaining task is to determine the parameter 
ranges over which reliable answers for amplitude statistics can be expected. 
Specifically, the level at which the variance of log amplitude saturates is usually 
regarded as the upper limit for useful amplitude calculations. The available 
experimental data on this subject are contradictory. Section 2.3 includes a discus- 
sion of these data and presents new data which may help to clarify the problem. 
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Little can be said about phase statistics, since experimental data is so 
sparse. We note, however, that second-order statistics such as the phase structure 
function are not appreciably different for the geometrical optics and the wave optics 
approach. In any case, it is thought that solutions for phase statistics using 

45 

Rytov's method will be valid for wider ranges than those for amplitude statistics. 

There are a number of well-known theoretical approaches which have been 

57 

neglected for various reasons. Hufnagel and Stanley's Mutual Coherence Func- 

2 

tion approach gives results which are in exact agreement with Tatarski ; in 
fact, their results may be considered a specialization of Tatarski' s more general 
results. The renormalization technique advocated by deWolf 52 * 65 * 6 ® has the in- 
teresting facet that it predicts a saturation of ir radiance fluctuations. In addition, 
the amplitude fluctuations are said to follow a Rayleigh or Rice distribution in 
place of the usually assumed log-normal behavior. However, the approach with 
its various approximations has not yet been subjected to a thorough critical re- 
view. In addition, the predicted saturation level is appreciably less than some 
of the measured values. 

In summary, it appears that Rytov theory, although it has certain limits 
which at this time are not well defined, is the most fruitful approach available 
to the propagation problem. Subsequent sections of this report will emphasize 
results obtained from Rytov theory. 

1.5 Theoretical Results 

This section is a summary of theoretical results which have been evaluated 
for various types of optical propagation problem. We restrict ourselves to 
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those results which have developed based on the Kolmogorov Theory of Turbulence. 
Therefore, we feel that the validity of contemporary optical propagation theory 
can be rigorously tested by comparing these results to the outcomes of properly 
designed experiments. 

The notational differences between authors is considerable, and we have at- 
tempted to standardize. This occasionally causes the constants in various equa- 
tions to differ slightly from the original published values. The classification of 
results is according to type of source (infinite plane wave, spherical wave, finite 
gaussian beam) and propagation path (homogeneous turbulence, non-homogeneous 
turbulence). Non-homogeneous means simply that C n 2 is a slowly varying func- 
tion of position along the path. The assumed spatial dependence will be given for 
each case. For each group of results, the listed order agrees chronologically 
with appearance in the published literature. The classification of results is as 
follows : 

(1) Infinite plane wave, homogeneous turbulence. 

(2) Spherical wave , homogeneous turbulence . 

(3) Finite gaussian beam, homogeneous turbulence . 

(4) Infinite plane wave, non-homogeneous turbulence. 

(5) Spherical wave, non-homogeneous turbulence . 

(6) Finite gaussian beam, non-homogeneous turbulence . 
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■ Common notation: 
Log amplitude 'C(x) 


LIST OF NOTATIONS 


= 1/2 In [u(x)/u Q (x)] , 

where u(x) is the field at x and u 0 (x) is the RMS field 
value at the same point. 

Optical carrier phase: />( x) . 

Refractive index structure constant as defined through 
the Kolmogorov "2/3 law": C n 2 . 

L and -£ : outer and inner scales of turbulence 

O O 

which define the limits of the Kolmogorov spectrum. 

Wave number of the radiation: k = 2n/K . 

Length of propagation path: L . 

Inner scale parameter: k m = (5.48/^ o ). 

Infinite plane wave: 

Log amplitude covariance: 

C* O) = <^(f) - ^(r))) (*( r + P) - (*<f + P)))y ■ 

where 7 and J> are position vectors perpendicular to the 
propagation path, and the assumptions of homogeneity 
and isotropy for 't have been made by writing (p) 

as (p). 

Log-amplitude variance: (0) = ^|^t(7) - (-£(7)^J ^ • 

Log-amplitude structure function: (p) = ^J/£(7) - £(7 + p)] 

= 2[cx(0 )-Qe(p)] 

Phase structure function: (p) = ^[/(r) - P(7 + p >]*>• 

Mutual coherence function: M(p) = (u(7)u* (7+p)}. 



■ iii 1 1 limn ii i I I 


Spherical wave: 

Log-amplitude covariance: (p) . 

Log-amplitude variance: C^ s (0) • 

Phase structure function: D^ s ( p) . 

Finite gaussian beam: 

Standard deviation of amplitude profile: a Q , 
i.e. , beam amplitude profile is proportional to 
exp [- £ 2 /2a 0 2 ] . 

On-axis log amplitude variance: Cy L ( 0 , 0). 

Log-amplitude variance as measured a distance R off the 
propagation axis: C^ L (0, R) . 

Mean value of log amplitude: (-C L ) . 

Phase structure function: ( s , S) , 

where s and S specify the position of the two measure- 
ment points with respect to the propagation axis. 

Theoretical Results 

a. Infinite Plane Wave, Homogeneous Turbulence . For/XL<<-i 0 (i.e., 
geometrical optics) , 

Q£(0) = 2.46C n 2 L 3 V /3 • (36) 

r- 

C^(P) = 277(0. 017) L 3 c n 2 J 0 (*p) k 4/3 ck, (37) 

•' 0 

D^(p) = 3.44k 2 LC n 2 p 2 I' 1/3 , with p «l 0 , (38) 
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and 




5 


D^(p) = 2.91k 2 LC n 2 p s/3 , with p » * B . 

For V\L » / & o (i.e. , wave optics), 

C^(0) = 0. 31 C n 2 k 7/6 L 11/6 , 


(39) 


(40) 


Clip) = 2rr 2 (0. 033) C n 2 k 2 L 
and 




k 2 L 


in — r— V /< 8/3 ck, 


(41) 


for p» -i , 

' O 7 

D^(p) = 1.72C n 2 k 2 L^e; 1 ' 3 P 2 , for p «l Q . (42) 


Also, for p >> -t 0 and for P - /XL , 

D^(p) = 2.91k 2 LC n 2 p 5/3 , 

and for ,o < /XL , 


(43) 


D^(p) = 1.46k 2 LC n 2 p 5/3 


(44) 


For p «% , 

D 0 (p) “ 1 • 7 2 C n 2 Zq 1/3 k 2 L p 2 . (45) 

2 

The above results, Eqs. 36-45, are due to Tatarski. The integrals in Eqs. 37 
and 41 have been evaluated numerically by him and are shown in his Figs. 8 
and 13. 

The following results (Eqs. 46-48) are due to Fried. 67 
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For /\L » / t Q , 


Qfc (p)/c.£(0) 


n= 0 


C) 

■3. 

■X: 


(k p 2 \ 2n / 

l^nWJ 


j /( 2n > ! 


/kp 2 \ 5 ' 6 

7. 53034 I , (46) 


Ci (0) = 0. 309 k 7/6 L 11/6 C n 2 , 


and 


D* (P)/Q(0)=2^- 1 + 7.53034 


'kp 2 \ s/6 ^->r (kp 2 \~| fkp 2 \ 2n 1 

4L J + ^[ a n +bn * 4L 4L / (2n)! 

n- 0 > 


(47) 


(48) 


The first ten values of a n> b n used in Eqs. 46 and 48 are listed in Table I and 
the results of Eqs. 46 and 48 are plotted in Fried's Figs. 1 and 3.®^ 

TABLE I 

Coefficients for Eqs. 46 and 48 


n 

a 

n 

b 

n 

0 

1.00000 

6.84209 

1 

-7.63889 x 10' 1 

5.27939 x 10~ 2 

2 

1.23778 x 10' 2 

-4.00354 x 10“ 3 

3 

-2.83086 x 10“ 3 

8.98375 x 10 ~ 4 

4 

1.08825 x IO - 3 

-3.09202 x 10 " 4 

5 

-5.34385 x 10' 4 

1.34605 x 10 ' 4 

6 

3.03065 x 10“ 4 

-6.80418 x 10 ~ 5 

7 

-1.89046 x 10' 4 

3.81508 x 10' 5 

8 

1.26184 x 10~ 4 

-2.30845 x 10 ' 5 

9 

-8.86012 x 10' 5 

1.48081 x 10 ' 5 
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The mutual coherence function for an infinite plane wave as measured on 


a plane perpendicular to the direction of propagation was evaluated by Hufnagel 
and Stanley- 57 as 


MO) = 



( 49 ) 


where the integration is over the propagation path. Therefore, for homogeneous 
turbulence , 

MO) ~ exp [- -j (2.91k 2 o s/3 C n 2 L )] • (50) 

b. Spherical Wave, Homogeneous Turbulence . From Tatarski, 2 for Y\L « -t o , 
C.£ (0) = 0. 246 C n 2 L 3 -O 7 ' 3 » (51) 

and , for V\L > > *t o , 

C^(0) = 0. 13C n 2 L u/6 k 7/6 . (52) 

From Fried, 58 for /XL>>-t o , we have 


Cf(p)/c^(0) = [^r] + e n [^] 17/6 }>^ 

n~ 0 ^ ^ 


ro 2 i 2n i r^ 2 i 5 ' 6 

[TO 698442 to * (53) 


C^ s (0) = 0. 124 k 7 ' 6 L 11/6 c n 2 , 


( 54 ) 
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and 



The first ten values of a n , b n , and c n are listed in Table II; the evaluation of 

£» q 

the right-hand sides of Eqs. 53 and 55 is shown in Fried's Figs. 1 and 2. 


TABLE II 

Coefficients for Eqs. 53 and 55 


n 

a 

n 

b 

n 

c 

n 

0 

1.00000 

6.84218 

2.53100 

1 

-2.59722 

-4.12854 

-0.389951 

2 

4.01006 

3.49205 

0.150525 

3 

-5.32819 

-3.17894 

-0.0790161 

4 

6.58359 

2.98285 

0.0485236 

5 

-7.79261 

-2.84418 

-0.0327814 

6 

8.96511 

2.73873 

0.0236145 

7 

-10.1076 

-2.65456 

-0.0178143 

8 

11.2249 

2.58505 

0.0139142 

9 

-12.3202 

-2.52616 

-0.0111669 


32 











c. Finite Gaussian Beam, Homogeneous Turbulence. From Kon and 


Tatarski 00 we have, for a collimated beam with /XL « -C , 


D^Cs, S) = 1.64 C 2 Z~ 1/3 Lk 2 s 2 


1 + 1 F 1 1/6, 2, 1.16 


S 2 L 2 k 2 
ka Q 2 


, (56) 


where 


and 


Pi " Pi 


s = (p, +P 2 )/ 2 • 

and p 2 are the position vectors of the observation points with respect to 
the beam axis. X F X is the confluent hypergeometric function, and the assumptions 
/XL « ■i o and near-field operation have been incorporated. In addition, Eq. 56 
has been written for the case where the observation points lie on a straight line 
passing through the beam axis. 

For the wave optics case (/XL >> -t 0 ) with L 2 k 2 / k 2 a Q 2 » 1, 


D* (s, S) 


0.617 k 3 a n C 2 k ~ 2/3 s 2 In 

0 n m 


(2. 7L 4 /< 4 

n 

^7 


[2 1 F 1 (1/6, 1, s 2 /a 0 2 )- jF, (1/6, 2, s 2 /a 0 2 )] • (57) 


From Fried 70 we have the following results for amplitude statistics. For 
collimated beam propagation , 


C^(0, 0)/c£ (0) = f. 


(£)■ 


(58) 
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where f- is evaluated for argument values ranging from 0.005 to 128.0 in Table HI. 


TABLE m 

Coefficients for Eq. 58 


( ka o 2/L ) 

f( ka o 2/L ) 

( ka o 2/L ) 


0.005 

0.9548 

0.60 

0.6237 

0.01 

0.9268 

0.70 

0.6590 

0.02 

0.8844 

0.80 

0.7000 

0.03 

0.8510 

0.90 

0.7442 

0.04 

0.8232 

1.0 

0.7901 

0.05 

0.7991 

2.0 

1.1947 

0.10 

0.7129 

4.0 

1.6120 

0.15 

0.6593 

8.0 

1.9239 

0.20 

0.6243 

16.0 

2.1356 

0.25 

0.6019 

32.0 

2.2712 

0.30 

0.5890 

64.0 

2.3708 

0.40 

0.5834 

128.0 

2.3774 

0.50 

0.5971 




For focused beam propagation , 


qf-co, 0)/cf (0) 



(59) 


where the function g is shown for discrete argument values in Table IV. 
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TABLE IV 


Coefficients for Eq. 59 


ka Q 2 /L 

g(ka 0 2 /L) 

ka Q 2 /L 

g(ka 0 2 /L) 

0.005 

0.9547 

0.60 

0.3691 

0.01 

0.9266 

0.70 

0.3402 

0.02 

0.8837 

0.80 

0.3157 

0.03 

0.8496 

0.90 

0.2948 

0.04 

0.8208 

1.0 

0.2766 

0.05 

0.7955 

2.0 

0.1723 

0.10 

0.7001 

4.0 

0.09827 

0.15 

0.6329 

8.0 

0.05212 

0.20 

0.5811 

16.0 

0.02606 

0.25 

0.5390 

32.0 

0.01234 

0.30 

0.5039 

64.0 

0.005583 

0.40 

0.4476 

128.0 

0.002304 

0.50 

0.4041 




We emphasize that Eqs. 58 and 59 hold only for measurement points on or near 
the axis of propagation. The results (Eqs. 58 and 59) are plotted in Fried's 


Fig. I. 70 


Results for the behavior of the log-amplitude variance off-axis are available 
71 

from Ho. These results are restricted to measurement points which lie in 
the near field of the transmitter aperture as well as to collimated beams. 

For YkL « -t o , 


qf-(0, R) 


c&(0) 



(60) 
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and 


C^( 0 ) 

= 0. 0108 T( 1/6 ) c n 2 K m 7/3 L 3 , 



/L 2 * n 2 \ 

c£{ o. o) 

= Qe(0) 



( 61 ) 


(62) 


For /\L » l o , 


C^(0, R) 

C-£(0) 

and 

0^(0, 0) 


c^(0) 


1-21- 


R' 


2\ / L 2 k 2 \ / L/< 


2 \- 5/6 


m \ I m 


a/ Ak 2 a 


2 „ 2 


2L 

ka, 2 


0. 0552 T( 1/6) C n 2 k 7/6 L 1V6 


C^(0) 


1 - 1.97 


- i o2 


(63) 


(64) 


(65) 


For clarity, we restate that R is the measurement point's off-axis radial distance 
and (0) is the log-amplitude variance for an infinite plane wave traversing 
the same path. Ho’s results are plotted in his Figs. 1 and 2. 71 

More general results are available from Ishimaru. 7 ^ These results are 
valid for collimated beams and a vanishingly small inner scale. The most useful 
results of this work take the form 

fa L 1 5/6 

q L (0, R) = tt 2 ( 0. 033 C n 2 )r(-5/6)k 7/6 L u/6 X 



2R 2 \ _ 

2a 0 2 [l+ (aL) 2 ]/ 


g(a-L) >• , 


(66) 
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where 


r 1 1 


i [l + (aL ) 2 x] [1 - 

xjl 5/6 

g(txL) - 1 Re < 

(1-x) 2 

aL 

— r dx 

•'o 1 



J 

6 J 

/ 1 + iaL 

\ S/6 ( 5 , 

17 \l 

- TT Re | 

V iaL 

) ' 2 F 1 \ 6 ’ 1 ’ 

6~ ; 


2 Fj is the hypergeometric function and a - K/ 2n a Q 2 . On the beam axis, the 
log -amplitude variance is 

</(0, 0) = 7r 2 (0.033 C n 2 ) r(-5/6)k 7 / 6 L 11/6 {1 ~ g(*L)}-(68) 


Eq. 68 is plotted as a function of aL in Ishimaru's Fig. 4, and it is shown there that 
C^ L (0, 0) increases initially at a rate between the rates for plane and spherical 
waves, eventually crosses the spherical -wave value, and thereafter grows at 
a lesser rate. 

The mean value for log amplitude has been evaluated by Gebhardt and Collins 73 
for the case where the detector is on or near the propagation axis and is in the 
near field of the transmitter. The analysis is based on prior work by Schmeltzer 74 
and is valid for focused and collimated-beam propagation. 

For collimated-beam propagation, it is seen that 


CO/ Qe (0) 



(69) 


where M was evaluated at discrete argument values as shown in Table V. 


37 



TABLE V 


Coefficients for Eq. 69 


ka Q 2 /L 

M (ka Q /L ) 

ka 0 2/ L 

M(ka 0 /L) 

0.005 

-2.796590 

0.60 

-2.071238 

0.01 

-2.730011 

0.70 

-2.089343 

0.02 

-2.632161 

0.80 

-2.108799 

0.03 

-2.558014 

0.90 

-2.128209 

0.04 

-2.497791 

1.0 

-2.146857 

0.05 

-2.447162 

2.0 

-2.273239 

0.10 

-2.276834 

4.0 

-2.365396 

0.15 

-2.180511 

8.0 

-2.417845 

0.20 

-2.122361 

16.0 

-2.445360 

0.25 

-2.087121 

32.0 

-2.459626 

0.30 

-2.066645 

64.0 

-2.468092 

0.40 

-2.052230 

128.0 

-2.473000 

0.50 

-2.057157 




For focused-beam propagation, 


U)/c£{0) = n (^) 


(70) 


This result is shown in Table VI, and both results (Eqs. 69 and 70) are plotted 
in Gebhardt and Collins' Fig. I.' 73 








TABLE VI 


0.005 


Coefficients for Eq. VO 


-2.796449 



N(ka Q 2 /L) 


-1.755759 


-2.729521 


-1.730162 


-2.630392 


-1.711608 


-2.554293 


-1.698131 


-2.491542 


-1.688406 


-2.437886 


-1.677912 


-2.246825 


-1.740934 


-2.123765 


-1.838889 


-2.035778 


-1.936671 


-1.969364 


-2.015523 


-1.917504 


-2.071282 


-1.842334 


128.0 


-2.107402 


-1.791461 


d. Infinite Plane Wave, Non-Homogeneous Turbulence. From Tatarski, 


we have, for geometrical optics (l^L << -£ 0 ) , 


Qi (0) - 7. 37^ o " 7/3 f C n 2 (x)x 2 dx, 

* n 


and, for wave optics (/kL >> -£ 0 ) , 


Ct (0) = 0.56k' 





where the integration variable x is the distance from the observation point to the 


part of the path being considered. Also , 


D*(/°) + tyC/o) = 2. 91 k 2 p 5/3 I C n 2 (x)dx 


J c n 2 (x; 

•'o 


(73) 


For 't << p << /\L, 


D, (p) = 1.46k 2 p 3 


' 3 j C’<* 

•'O 


) dx 


(74) 


For p - V\L >> -i , 


D*(p) = 2.91k 2 p 5 / 3 C n 2 (x)dx 

‘o 


f c - 

J 0 


(75) 


Two models for the turbulence are used for calculations. For 


C n 2 ( X ) " C n0 eX P [" X/X 0 ] ' 


the wave optics solutions become 


C£ (0) - 0. 53C n 2 0 exp [- x./x 0 ] 


and 


% 


2.91 k 2 p 5/3 C_ 2 „ x 


n0 "0 ’ 


(72A) 

(75A) 
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and, for 


they become 


and 


C 2 (x) 

n v ' 


"n 0 



( x/x o) 2 J 


C^(0) = 3.4C„ 2 k 7/6 x Q 11/6 

D*(P) = 4- 57 C n 2 k 2 p 5/3 x Q . 


(72B) 


(75B) 


The above results (Eqs. 71-75B) are discussed in the context of the stellar 
"seeing" problem in Chapter 8 of Tatarski.2 Fried^ 6 has investigated the same 
problem using the more realistic turbulence model 

C n 2 (h) = 4.2 x 10‘ 14 h" 1/3 exp [-h/h 0 ] , (76) 

where h is the altitude measured from sea-level and h 0 is the scale height of 
the atmosphere (3200 meters). The results are 


q> (0) = 0.730 cos 




(77) 


V - ' 1 f 


Qfc(p)/Ct(0) = 


k p 2 


4h Q sec^ 


\. ( kp 2 y /2 / k/p 2 tvx 

J c n \4h 0 seciS 1 j \4h Q seed 


(4h7^r. <™> 
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and 


(V)/'C^(0) 


2 <- 1 + 6. 27602 


ko 2 

4h 0 sec 0 


5/6 


OD 

L 

n= 0 


+ b„ 


k p 7 


4h Q sec 9 


+ c 


kp 5 


4h Q sec 6 


3/2 


k p 2 


, 5 - 2 


4h Q sec 9 , 


(kP 2 /4h 0 


[ ( 2n ) ! ] 



(79) 


where 9 is the zenith angle of the source and the coefficients a , b , c ,and d 

° n n n n 

are listed in Table VII, The results of Eqs, 78 and 79 are plotted in Fried's 
Figs. 2 and 4. 7 ® 

TABLE VII 

Coefficents for Eqs. 78 and 79 


n 

a 

n 

b 

n 

c 

n 

d 

n 

0 

1.00000 

6.22008 

-7.85622 x 10" 1 

1.45145 x 10~i 

1 

-5.55556 x 10 _1 

-1.79179 x 10- 1 

2.28026 x 10“ 2 

-1.73367 x 10 - 3 

2 

3.74486 x 10-1 

5.05797 x 10"2 

-4.36773 x 10' 3 

1.67119 x lO’ 4 

3 

-3.13723 x 10-1 

-2.24463 x 10“ 2 

1.45802 x 10" 3 

-3.32886 x 10“ 5 

4 

2.79596 x 10" 1 

1.23097 x lO- 2 

-6.40001 x 10" 4 

9.68214 x 10" 6 

5 

-2,56692 x 10 1 

-7.63945 x 10" 3 

3.30962 x 10“ 4 

-3.55757 x 10~ e 

6 

2.39822 x 10" 1 

5.14252 x lO" 3 

-1.90915 x 10" 4 

1.53244 x 10" 6 

7 

-2.26661 x 10 - 1 

-3.66677 x 10" 3 

1.19079 x 10" 4 

-7.40745 x lO - 7 

8 

2.15982 x 10 _1 

2.72914 x lO" 3 

-7.87611 x 10" 5 

3.90776 x 10" 7 

9 

-2.07068 x 10" 1 

-2.09991 x lO" 3 

5.45286 xlO" 5 

-2.20795 x lO" 7 
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e. Spherical Wave, Non- Homogeneous Turbulence. From Fried 77 we have 


the following general result: 


. f L , 

x\ 5/3 


D^(P)+D^(P) = 2.91 k 2 p 5 ' 3 I C n 2 (x)^ 

•'o 

r) dx. 

( 80 ) 


where, in this case, x is the distance from the source to a point on the propaga- 

70 

tion path and L is the total path length. Also, from Fried we have 


q?(0) 


0. 727 sec 11 ' 6 9 



(81) 


where the turbulent profile of Eq. 76 has been used and it is assumed that the 
measurement point is at a much greater altitude than h Q (3200 meters). 

70 

f. Finite Gaussian Beam, Non- Homogeneous Turbulence. From Fried, 


c^(0, 0 )/c| ( 0 ) = 7"( a „ +b n fl+ c n n 3/2 + d n n 5/2 )x 

n= 0 


/ fi 2n \ 

7 — - 5. 90358 fi 5/6 > 

\[ 2 n ] ! 1 


(82) 


where 


n 


ka o 2 


h Q sec . 
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and the coefficients are listed in Table VIII. Again, the turbulence model of Eq. 76 
has been used and assumption of measurement ranges much greater than h Q has 
been applied. Eq. 82 is plotted in Fried's Fig. I. 78 


TABLE Vin 
Coefficients for Eq. 82 


n 

a 

n 

b 

n 

c 

n 

d 

n 

0 

1.00000 

6.22009 

-1.04436 

4.82369 x 10" 1 

1 

-5.55555 x 10" i 

-5.37538 x 10" 1 

1.32617 x 10" 1 

-4.53727 x 10" 2 

2 

3.74486 x 10" 1 

2.52899 x 10" 1 

-5.23919 x 10" 2 

1.30301 x 10" 2 

3 

-3.13723 x 10" 1 

-1.57124 x 10“ 1 

2.84200 x 10" 2 

-5.51539 x 10“ 3 

4 

2.79596 x 10" 1 

1.10788 x 10" 1 

-1.79886 x 10" 2 

2.85744 x 10" 3 

5 

-2.56692 x 10" 1 

-8.40340 x 10" 1 

1.24807 x 10" 2 

-1.67697 x 10" 3 

6 

2.39822 x 10" 1 

6.68528 x 10" 2 

-9.20387 x 10" 3 

1.07123 x 10" 3 

7 

-2.26661 x 10“ 1 

-5.50016 x 10" 2 

7.08917 x 10" 3 

-7.27633 x 10 " 4 

8 

2.15983 x 10" 1 

4.63955 x 10" 2 

-5.64133 x 10" 3 

5.17810 x 10" 4 

9 

-2.07068 x 10" 1 

-3.98984 x 10" 2 

4.60446 x 10“ 3 

-3.82207 x 10" 4 


The preceding equations (Eqs. 36-82) are, to the author's knowledge, a 

complete listing of numerically useful results now available based on Kolmogorov 

turbulence theory. At the time of this writing, it appears that most theoretical 

efforts are being directed toward the finite gaussian beam propagation problem. 

This approach is particularly interesting since it is compatible with experiments 

utilizing laser sources. Also, the finite beam solutions must reduce in limiting 

cases to the spherical -wave and infinite -plane -wave results. The work of the 

79“ 8 1 

Japanese group of Kinoshita, Asakura, and Suzuki has been particularly 
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noteworthy. Using the Rytov approximation and a gaussian model for the covariance 
function of the refractive index field, this group has evaluated various second- 
order statistics for — 

(1) Log amplitude. 

(2) Phase. 

(3) Joint log-amplitude-phase fluctuations throughout the beam. 

The results pertain to homogeneous turbulence and focused, collimated, or 
uncollimated propagation. We do not list these results, however, since their 
accuracy is questionable owing to the gaussian turbulence model employed. 



CHAPTER 2 


EXPERIMENTAL RESULTS 


2. 1 Experimental Program 

a. Introductory Remarks . The experimental side of the optical propagation 
problem has been investigated in one form or another for many years. Historically, 
astronomers who have studied the brightness and position fluctuations of stellar 
images must be regarded as the first experimentalists. However, these workers 
often failed to recognize the problem in its true statistical form; their meas- 
urements are difficult to interpret and, therefore, equally difficult to relate to 

2 

contemporary theories. The pioneering theoretical treatise by Tatarski also 
contained a substantial amount of experimental work, but these results have been 
criticized in recent years as regards the type of optical source and the data 
analysis procedures used. In general, the experimental situation is very con- 
fusing; only a few published results are widely accepted by workers in the 
field. 

It should be clear that optical propagation measurements are useful only if 
they are taken in such a manner as to permit direct comparisons with theory. 
Otherwise, the measurements must be regarded as pertinent only for the time, 
location, and meteorological conditions under which they were taken. In principle, 
the experimenter should have knowledge of the statistical state of the turbulence 
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along the entire propagation path. Then a wave of the proper type (e. g. , spherical, 
infinite plane) is launched into the medium, and the received fields are measured 
and statistically analyzed at the receiver. Practical limitations dictate that the 
number of meteorological stations be small, and one therefore usually chooses 
paths which are nearly uniform topologically. The measurements at a single 
station or the average of measurements at a few stations are used to characterize 
the turbulent state of the atmosphere along the entire path. 

There are two primary reasons for seeking a theory which accurately pre- 
dicts the effects of turbulence on optical waves. In the first case, the efficient 
implementation of optical systems for informational purposes requires that the 
designer be able to predict accurately the optical carrier fluctuations induced by 
the known (or assumed) turbulence. Therefore, for known turbulence models, 
systems can be optimized to reduce the effects of turbulence. In the second case, 
the optical fluctuations are used to infer the state of the turbulence. The meas- 
ured field quantities are substituted into the theory, and a unique solution for the 
turbulence which caused the fluctuations is sought. Either approach requires the 
existence of an accurate theoretical treatment of the propagation problem. 

In the work that follows, we emphasize the first type of problem, the optical 
communication problem. The experiments are designed to measure those param- 
eters which control the quality of a communication channel. In all cases, the 
measurements are carried out in a manner which permits direct comparison with 
theory. The results to be presented hold for near-ground horizontal paths. 

b. Test Facilities . All of the experiments to be described were performed 
on a test range operated by Goddard Space Flight Center about 5 miles east 


48 



of Greenbelt, Md. The range is located in an agricultural research area and con- 
sists mostly of live-stock graze-lands. The optical path lies approximately in 
an east-west geometry and has winds which are usually out of the northwest. 

At the west end of the path, a permanent structure houses all transmitting, re- 
ceiving, and recording equipment. Two large mirrors are located on the range 
and serve to reflect transmitted beams back to the transmitter station. The 
transmitter and receiver systems are separated by about 1 meter, so that the 
optical path to and from the mirror has overlap only in a small region near the 
"folding" mirror. A profile of the optical path is shown in Fig. 1. The average 
beam height about ground level is 4.2 meters with a standard deviation of 
1.3 meters. 

The transmitter-receiver building is a 20' by 20' structure (Fig. 2) which 
is partitioned into two 10' by 20' sections. During experiments, one of these 
sections is well enough ventilated so that conditions remain near ambient. The 
other section contains all the apparatus requiring a stable environment and 
is heated and air-conditioned. The transmitter and receiver systems are mounted 
on a one-ton granite slab (Fig. 3). This table has three supports, which are iso- 
lated from the building floor and rest on the eight-inch-deep concrete pad beneath 
the building. The granite table can be "floated" on pneumatic supports (com- 
pressed nitrogen) for additional vibration isolation. In this condition, the system 
rejects input excitations of frequency 1 Hz or greater. 

The first folding mirror is 50.8 cm in diameter and is located at a range of 
583 meters. The mirror is in a two -axis adjustable mount (Fig. 4) and is attached 
to a one-ton concrete block. The second mirror (Fig. 5) is 30.5 cm in diameter 
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Figure 2. Transmitter-receiver station. The build- 
ing is 20' by 20' and sits on an 8 -inch-deep concrete 
pad. Transmitted and received beams both go through 
the open window shown above. 



Figure 3. Stable transmitter-receiver platform . The 
compressed nitrogen tanks in the background can be 
used to "float" the table for additional vibration 
isolation. 


I 
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and is also in an adjustable mount and located at a range of 200 meters. Both 
mirrors are optical -quality flats and are protected by shelters when not in use. 

During data runs, the signals of interest were recorded on magnetic tape. 

The recorder was an Ampex FR-600, which is a 7-channel machine having 
either direct- or fm-recording capability on each track. Binary coded decimal 
time signals were also recorded, so that efficient data processing could be imple- 
mented. In all cases the recorder was run at 7-1/2 inches per second, which 
resulted in an fm bandwidth of DC to 2.5 kHz. Figure 6 shows the tape recorder 
and associated signal electronics. 



Figure 6. Magnetic tape recorder and supporting 
electronics. This system was located in the heated/ 
air-conditioned portion of the building and was con- 
nected by about 25 feet of cable to the optical 
detectors. 
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Several data-processing techniques were used, depending on the statistic of 
interest. We discuss here only the digital computer approach and note that the 
other techniques will be discussed as they occur in the sections to follow. Two 
computers were used during the course of the program. The first was a Raytheon 
520, which is a relatively small machine (8K memory). This machine had input 
analog-to-digital converters, with sample rates variable up to 10 4 per second. 

The playback of the tape recorder was sampled by these units 10 3 times per 
second and stored in the memory in correct time sequence. Data analysis was 
initiated manually by a switch and terminated when the designated number of 
samples had been stored. The program analyzed both the amplitude statistics 
of the signal (i.e., probability density, cumulative density, moments) and the time 
statistics (autocorrelation, power spectral density). After this procedure, each 
stored data point was normalized by the mean value of the original signal ; then 
the natural logarithm was taken to produce the variable log amplitude. Identical 
amplitude and time analysis was then performed for log amplitude. 

In order to process longer time records it was necessary to use a larger 
machine. The IBM 360/91 was therefore programmed to do the same analysis 
as outlined above, on longer data records. In this case, the analog data tape had 
to first be converted into a digital data tape. The digitization was carried out at 
a 1-kHz rate and with sufficient bits per sample point so that quantization noise 
was insignificant compared to other noise sources. 

During all experiments the average meteorological conditions (temperature, 
pressure, relative humidity, wind magnitude, and direction) were monitored at 
the transceiver station. In addition, a pair of high-speed microthermal probes 
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were mounted near the west end of the path at a height of 4.15 meters in order 
to make direct measurements of the strength and spectrum of the temperature 
fluctuations. Figure 7 shows the pair of thermal sensors mounted on a wind- 
vane -type structure which was used to keep the sensor elements perpendicular 
to the wind vector. The center probe is a velocity sensor which was used for 
accurate wind speed measurements. The thermal sensor elements were 2.0- 
micron-diameter strands of platinum with a length of 1.6 mm. The system had 
a resolution of 0.01°C and a bandwidth of 1 kHz. The processing electronics for 
these sensors is shown in Fig. 8. 



Figure 7. Micrometeorological meas- 
urement system. The center probe 
measures wind speed, while the two 
outside probes monitor microthermal 
fluctuations. The system was mounted 
4. 2 meters above ground, which was 
the average height of the propagation 
path. 


55 




Figure 8. Processing electronics for 
micrometeorological system. The 
unit at the top of the rack is the ve- 
locity meter. The remainder of the sys- 
tem processed temperature and tem- 
perature -diffe rence fluctuations . 


We note that this experimental program was approximately a 2 -year effort 
starting in Dec. 1967. Not all of the instrumentation described in this section 
was available at all times during the program. For example, the microthermal 
sensor system did not become operational until near the end of the program. 

Also, the shorter range capability provided by the 30. 5 -cm mirror was not avail- 
able at the time of the dual -wavelength experiment (Section 2.2). 
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2.2 Wavelength Dependence of Laser Beam Scintillation* 


a. Introduction . It is shown in Appendix B that the variance of log amplitude 
calculated through Rytov theory has a wavelength dependence k (3/2 ” p/2) , where 
p describes the spatial behavior of the refractive index fluctuations through 

( [n(x + r)-n(x)] 2 ) ~ r p . (83) 

Using the Obukhov-Kolmogorov turbulence model, we obtain 

P = 2/3 . 

Therefore the strength of optical fluctuations (as denoted by the log-ampli- 
tude variance) can be expected to vary as the 7/6 power of the wave number. The 
experiment described in this section was designed to test this theoretical 
prediction. 

b. Experimental Technique . Because of the short-term stationarity and 
lack of statistical homogeneity of the atmosphere, the experimenter must trans- 
mit simultaneously through the same portion of the atmosphere at each wavelength 
of interest. The experiment is shown schematically in Fig. 9. A He-Ne laser at 
0.632 micron and a C0 2 laser at 10.6 microns were aligned parallel, with an 
offset of 12 cm, and were operated over a folded horizontal 1.17-km path. The 
He-Ne output was unmodulated, and no external optics were used; thus, the 
assumptions of a spherical -wave transmitter were satisfied. The C0 2 output 

*A substantial portion of the results of this section was presented to the Spring 1968 meeting of 
the Optical Society of America and subsequently published 82 in the Journal of the Optical Society 
of America. 
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APERTURE STOP 



Figure 9. Experiment schematic diagram, dual -wavelength experiment. 


was modulated at a 1-kHz rate to minimize the background and to enhance the 
detector operation. The 10.6-micron beam was transmitted through a 13-mm- 
aperture telescope which, because of the longer wavelength, also approximated 
a spherical -wave source. A small portion of the C0 2 beam was diverted and de- 
tected by a thermopile to monitor any fluctuations in output power. 

A mirror in the field reflected the two beams, and they were then detected 
at the transmitter station by a photomultiplier and a gold-doped germanium de- 
tector. The output of the photomultiplier tube was recorded directly on an FM 
tape recorder. The 10.6-micron detected signal was demodulated by means of a 
standard envelope detection scheme and then recorded. 
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Table IX lists the important parameters of the experiment. In order to 
avoid the effects of aperture-averaging, each detector aperture was made to ap- 
proximate a point detector. In addition, the photomultiplier tube aperture was 
adjusted so as to make the ratio of the apertures equal to the square root of the 
ratio of the wavelengths. This ensured that aperture-averaging effects were 
insignificant. 


TABLE IX 

System Parameters, Dual -Wavelength Experiment 


Parameter 

O 

O 

to 

He-Ne 

Output power 

2 W 

1 mW 

Beam divergence 

0.4 mr 

0.7 mr 

Beam diameter 

13 mm 

1.4 mm 

Detector aperture 

2.5 cm 

0.62 cm 

Detector field of view 

2° 

5° 

Detection signal bandwidth 

220 Hz 

2.5 kHz 

Chopper frequency 

1 kHz 


Range 

1.2 km 

1.2 km 

Average beam height 

4m 

4m 

Folding mirror diameter 

50.8 cm 

50.8 cm 

Length of data run 

10 sec 

10 sec 


The recorded scintillation data were replayed through an analog-to-digital 
converter into a Raytheon 520 digital computer that analyzed the statistics of the 
scintillation. The computer was programmed to analyze both irradiance and 
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log-amplitude statistics. Available statistics and functional relationships are 
given in Table X. The analog scintillation signal was sampled at a 1-kHz rate 
for 10 seconds, with each sample assigned to one of 500 discrete levels available 
in the computer. The selection of the time-record length for a single data run 
is a tradeoff between the statistical reliability of a measurement and the statisti- 
cally time-varying atmosphere. 


TABLE X 

Computer Analysis, Dual -Wavelength Experiment 


Irradiance Statistics 

Log-Amplitude Statistics 

Average value 

Average value 

Variance 

Variance 

Median 

Median 

3rd central moment 

3rd central moment 

4th central moment 

4th central moment 

Probability den sit} 7 function 

Probability density function 

Power spectral density 
Autocorrelation function 

Cumulative distribution function 


c. Results . Fig. 10 shows results of measurements taken over a 2 -hour 
nighttime period in December. Plotted along the Y-axis is the parameter , de- 
fined as the ratio of log-amplitude variances at 0.632 and 10.6 microns. The 

2 

X-axis is the log-amplitude variance for the He-Ne beam (o-^ 0 632 ) and the 
refractive index structure constant C n 2 , which in this case was determined from 
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cl X 10 1S (METERS) 

Figure 10. Ratio of log-amplitude variances for He-Ne and C0 2 wavelengths. 

the equations for horizontal spherical -wave propagation. The value 26.7 is the 
theoretical ratio obtained from the 7/6 power predictions. 

The data would seem to indicate that the scintillation statistics depart from 
the 7/6 power law in the region to the left. The difficulty here is that the region to 
the left corresponds to very low turbulence. The C0 2 laser output unfortunately 
had an appreciable amount of 120-Hz power supply ripple. In this region, the 
dominant noise on the 10.6-micron signal was caused not by the atmosphere but 
by the 120-Hz ripple. The output of the C0 2 unit was recorded on tape; an 
evaluation of the magnitude of the ripple showed that, in the 0.04 region, the 
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internal laser modulation became the same order of magnitude as the atmosphere - 
induced modulation. Therefore, the data points in the region below 0.07 are of 
questionable value. Some compensation for the C0 2 internal noise can be ac- 
complished if one makes approximate assumptions. 

The intensity at the detector aperture can be considered as the product of 
two random variables. 


I(t) = I L h(t) . 

The first, I L , is the irradiance fluctuation caused by internal laser ripple; the 
second, h(t), is the time-varying atmospheric transmission. If both variables 
are assumed log-normal with respective variances o^ 2 and o- 2 2 , their product is 
also log-normal with a variance equal to o + o- 2 2 . 83 Several experiments 70 * 84 
have indicated that h( t ) is log-normal in the visible spectrum, and inspection of 
cumulative density of log amplitude for the transmitter noise (Fig. 11) has in- 
dicated normality as a reasonable approximation. 

The data are shown replotted in Fig. 12 with the transmitter noise subtracted 
from the C0 2 statistics; the measurements obtained during low- scintillation con- 
ditions are neglected. The average value of the 39 independent data points shown 
is 26.8, and their variance is 10.8. The average value of Jj is in close agreement 
with the predicted value of 26.7. 

The statistical reliability of the measurements has been analyzed to de- 
termine the magnitude of data scatter to be expected under ideal conditions. 

Consider two target populations x and y where the random variable x and y 
are independent gaussian random variables. (Note that x and y represent the 
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LOG AMPLITUDE 



0.10.5 2 5 10 20 40 60 80 90 95 99 99.8 99.99 

CUMULATIVE PROBABILITY 

Figure 11. C0 2 transmitter noise; log-amplitude 
cumulative probability. 


« 26.7 
25 


*=26.8 
VAR[*]= 10.8 


Cnx 10 15 (METERS) 

Figure 12. Ratio of log-amplitude variances for He-Ne and CO, wavelengths 

2 ^ 

(transmitter noise removed, and all data o-^ Q 632 less than 0,07 neglected). 




variable log amplitude at 0.632 and 10.6 microns; the assumption that they are 


independent is thus a very good approximation, as shown later in Fig. 32.) Let 


and 


VAR(x) = cr x 2 
VAR(y) = a 2 . 


(84) 


Select a random sample of size N from x and form the sample variance 


S 


2 

X 


1 

N~ 1 




(85) 


Since the x i are normally distributed, then the random variable 


(N- 1) S x 2 


has the probability density function 

(x 2 ) N+V2 

p (* 2 ) = ex P [- X 2 /2 ], for x 2 >0, (86) 

2 N ' 1/2 r( il T ± j 


which is the chi-square density with N - 1 degrees of freedom. So, in this case, 

2 - 2 
X " x N -i • 
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Similarly, one can repeat the experiment with the y population and form the 
sample variance 


V = F-r 2j y ‘“ 5 ' 2 ' (87) 

i = l 

where the random sample is of size M. Note that the quantities N and M refer to 
the number of independent samples and not the total number of samples. 

As before, (M - 1) S y 2 /cr y 2 is a chi-square random variable withM - 1 degrees 
of freedom. Therefore , 

_ (M-l)S y 2 

X M -i “ ~ ~2 ' (88) 

y 

We now consider the normalized ratio of the sample variances, i.e. , 

S 2 /<r 2 

XX 

s 2 /a 2 = F n-i.m-i • 

y y 

This quantity is a function of a random variable and is therefore itself a random 
variable. It is well known 85 that F n _j M-1 has the density 


p( f n-i.m-i) 


/N+M \/N-l) rl/2 , . 

\ 2 V \ M~ 1/ ( F n-i. m-i) 

~ i\ jm - i\ r fN - 1\ 

l 2 ) r \ 2 ) L 1 + \M - 1/ F N-1,M-1 


N-3/2 


N+M-2/2 


for F — 0 ■ (89) 


So, in order to construct a 90% confidence interval for a particular set of meas- 
urements, we write the probability statement 

P [ a<F n-i,m-i < B I = 1_a = - 90 ’ (90) 
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mu ii i mu i ii 1 1 1 ii 


where B can be determined directly from tables, once N and M are known. In order 
to find A , we use the relation 


l N- 1 ,M-1 : l-a/2 


M- 1 , N- 1 : a/2 


(91) 


Therefore, 


L M“ 1 , N” 1 : a/ 2 


S 2 a 2 
— — — < F 

s 2 g 2 1 N-l,M-l:a/2 
y x 


o 2 p 

L y r N*l , M _ 1 :a/2 


s 2 f, 


x 1 M-l ,N-1 :a '2 


1 - a = .90 . (92) 


For the sake of interest consider the following situation: Assume that a time 
record of length T has been analyzed for both the x and y target populations and 
that the evaluated ratio of the respective variances is 27, i.e. , 



27 . 


hi order to determine the degrees of freedom for S x 2 and S y 2 , we inspect a typical 
autocorrelation function for each process (Fig. 14) and use a linear fit to the 
autocorrelation functions as well as the results of Appendix A ; the number of 
independent samples N and M is then evaluated as 


N = aT % 120(10) = 1200, for S x 2 , 

and (93) 

N = aT * 28.6(10) = 286, for S y 2 . 
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Therefore, from F-distribution tables (with linear interpolation) we get 


Therefore, 



a 

< —7 < 27(1.18) 

cr * 


y 


.90 . 


a* 

P 22.9 < < 31.8 

cr * 


.90 . 


( 94 ) 


*In a strict sense, the equivalent bandwidth used for finding N should be the bandwidth of the log- 
irradiance signal and not the bandwidth of the irradiance fluctuations. However, the data reduction 
system did not permit an evaluation of the sprectrum of log irradiance; therefore, the equivalent 
bandwidth of the irradiance signal was used. It is easy to show that the bandwidth of the irradi- 
ance signal is not seriously affected by the logarithmic operation under weak scintillation condi- 
tions (which is the situation at 10.6 microns, for most paths). Consider a system with input 
x = v Q = v(t) and output y = A log x, where 

Vq = constant (dc level), 

v(t) = random time function with mean value equal to zero, 

and 

A = constant (gain of the log amplifier). 


Power spectral density (psd) analysis of the input x is performed on v(t) only, since Vq merely 
results in a delta function at the origin in the frequency domain. The output of the system is 


y - A log [v 0 + v(t)] - A log (v Q ) 1+-^ - A log (v Q ) + A log 


1 + 


v(t) 


Therefore, psd analysis of the output needs to be performed only on the second term, since A log 
[v Q ] is a constant. However, 


log (1 t x) % x, for x << 1 ; 


therefore, 


log 



v(t) l 

v 0 J 


v(t) 


for v(t) « Vq (i.e., weak scintillation). Thus, psd analysis of y is equivalent to psd analysis of 
A v(t)/vg; this is the same as psd analysis of the input x, since the constant factors A, Vq are 
removed by dividing through by the mean-square value. 
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The length of the confidence interval for the other data runs can be expected to 
be about the same , since T = 10 seconds for all cases. In addition, B * 
is determined chiefly by the component of wind velocity perpendicular to the 
propagation path, and this did not change significantly during the experiment. 

So we conclude that the observed data scatter is of the same order of magnitude 
as could be expected from statistical variation. 

The normalized power spectral density for both wavelengths for a particular 
data run is shown in Fig. 13. These functions were evaluated by means of a 
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Figure 13. Dual -wavelength normalized power spectral density (resolution band- 
width = 5 Hz; each curve normalized to the same area; log-amplitude variances 
for these data are cr£ 0 632 = 0.200 and o-£ 1Q g = 0.00693). 
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Fourier transformation of their respective autocorrelation functions. Both curves 
are normalized by their intensity variances. Thus, if each wavelength has equal 
amounts of scintillation, this plot would indicate spectral distribution of the 
scintillation. The normalization scheme used requires the height of the C0 2 plot 
to be increased by a factor of 16. Therefore, it should be clear that, at any fre- 
quency, there is far more scintillation in the visible spectrum than at 10 microns. 

The resolution bandwidth is 5 Hz, and the wind velocity perpendicular to the propaga- 
tion path is approximately 1.5 meters per second. The power supply ripple on 
the 10-micron beam is obvious in this plot. 

Fig. 14 presents the same data in the time domain in terms of the normalized 
autocorrelation function R x (r)/R x (0) . Again, the internal modulation of the 
C0 2 laser appears as peaks at 8- and 16-millisecond intervals. If one fits an 
exponential function of the form e~ aT to each of these curves, the correlation 
time constants (the alphas) have a ratio of 4.4, which is close to the ratio of the 
square roots of the wavelengths. This is in agreement with propagation theory, 
as can be seen from the following argument. Tatarski has shown that the covariance 
of log amplitude (p) falls off to zero for separations on the order of ikz. Also, 
it is shown analytically in this report and experimentally by others 2 that the ir- 
radiance covariance function Cj (p) has a similar characteristic. So, (p)^0 
at P = Assuming that the relationship between the spatial dependence and 

the time dependence of the covariance functions is predicted by the concept of 
"frozen-in turbulence," then p - V n r, where V n is the component of wind velocity 
perpendicular to the propagation path; therefore, 
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Figure 14. Normalized autocorrelation function for dual -wavelength scintillation; 
fitting the data with the exponential e -aT gives a (C0 2 ) = 71.4, a(He-Ne) = 312. 

So, if we inspect C r (r) for two different wavelengths, we expect the ratio of 
the correlation times to be proportional to 2 • 

A typical portion of the analog scintillation signal is shown in Fig. 15. The 
computer evaluated the variances as 0.200 and 0.00668, which gives a 0 value of 
30. The chart recorder had a 110-Hz, 3-db bandwidth, which tends to smooth 
some of the high-frequency He-Ne scintillation (very little C0 2 scintillation is 
in the 110-Hz region). The important parameter is the ratio of the maximum 
and minimum signal excursions from the average value. For the He-Ne scintilla- 
tion, excursions of several octaves are common. In both of these traces, the 
background levels were not significantly different from the zero level. 
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Figure 15. Typical analog scintillation signal; log-amplitude variances are 
0.200 for He-Ne and 0.00668 for C0 2 . 

Fig. 16 is a typical cumulative probability plot of log amplitude at 10.6 
microns on Gaussian probability paper. The abscissa is F(i) , where 

f' c 

F(-C) - f(u)du , 


and f(u) is the probability density function for log amplitude; the ordinate is 
log amplitude. The abscissa scale is constructed so that a normal random 
variable would plot as a straight line. Thus it is clear that the statistics of 
scintillation at 10.6 microns satisfy this test for log-normality. 
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LOG AMPLITUDE 



CUMULATIVE PROBABILITY 

Figure 16. Cumulative probability of log amplitude at 10.6 microns; 
data plotted on a Gaussian abscissa scale. 


d. Conclusions . In summary, the following conclusions can be made: 

(1) The average value for the log-amplitude variance ratio at 10. 6 and 0. 632 
microns is very close to the predicted value of 26. 7. 

(2) The C O z scintillation has a bandwidth that is several times less than for 
the He-Ne and a correlation time that is approximately 4.4 times greater 
than for the He-Ne. 

(3) The probability distribution of scintillation at 10. 6 microns is approxi- 
mately linear when plotted on log-normal probability paper and thereby 
satisfies this test for log-normality. 


72 



2.3 Measurement of Log-Amplitude Variance and the Saturation Level of Irradiance 
Fluctuations* 

a. Measurement of Log-Amplitude Variance . Theoretical expressions for 

2 

the variance of log amplitude (oy ) have been derived for several types of propa- 
gation paths. 2,70,78 Numerous experiments have been devised to check the 
theory and have resulted in several techniques for the experimental evaluation 
of . The first part of this section is concerned with numerical differences 
in <?£ which we have noted in our work when two different techniques were applied 
to the same raw data. 

It has been shown 86 that for log-normal scintillation, irradiance and log- 
amplitude variances are related through 

a i 

7T = ex P ] ~ 1 . (95) 

where c ^ 2 = irradiance variance. Since irradiance fluctuations are easily meas- 
ured, several experimenters 2 ’ 84, 87 have used Eq, 95 to evaluate oy 2 . Another 
technique is to simply normalize a group of quantized irradiance samples (I lt 
I 2 , • ■ • I n ) by I 0 and form ('C 1 , • • • l n ), where 

, 1 , 

-u . - "9 In y- , for i - 1,2 n . 

A o 


*A portion of the work reported in this section was published in the Journal of the Optical Society 
of America.®® 
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A histogram (i.e. , probability density function) can then be constructed and the 
moments of log amplitude computed. 

Both of these techniques were used on scintillation measurements with a 
He-Ne laser (0.632 microns) over a 1.17-km horizontal path. The laser output 
was 1.4 mm between l/e 2 points in the irradiance profile and had a divergence 
of 0.7 mr. At the receiver, the fluctuating signal was detected by a photomulti- 
plier through a 3-mm aperture. The results are shown in Fig. 17, where o- T 2 is 
the log -amplitude variance based on irradiance measurements and Eq. 95, and 
c 7 £ 2 is the log-amplitude variance evaluated from the moments of log amplitude. 



Figure 17. Comparison of log-amplitude variance: cr T 2 calculated from irra- 
diance moments [Eq. (95)] and cr E 2 calculated from log-amplitude moments. 


74 



Ideally, either technique should give the same result and all points should 
lie on the unity-slope straight line in Fig. 17. However, as the scintillation be- 
comes stronger, the difference (o- E 2 - cr x 2 ) becomes significantly greater. There- 
fore, the discussion must center on which of the two techniques for <?£ - a T 2 
evaluation is more valid. 

Objections to the use of moments to evaluate cr E 2 are based primarily on the 
dynamic range and frequency response limitations of the instrumentation. In 
this experiment the average photomultiplier current was approximately a micro- 
ampere and the signal -to-noise ratio was sufficiently high so that dynamic range 
constraints were not excessive.* Power spectral density analysis of the data 
indicates that the frequency response of the system (2.5 kHz) was also not a 
problem. 

Since o" T 2 is based on the Gaussian behavior of log amplitude (through Eq. 95), 
it would seem reasonable to try to correlate (<y£ - cr*) with deviations from 
normality. Fig. 18 plots percent deviation versus skewness S, where S = m 3 / 

( m 2 ) 3 2 and nr is the ith moment of log amplitude about the mean. For a sym- 
metrical distribution such as the normal, S should vanish. Fig. 18 is shown with 
a least- square linear fit. 

The data indicates that for scintillation with | s| >1, the error in calculating 
log-amplitude variance through Eq. 95 becomes significant. If Fig. 17 is extrap- 
olated to stronger turbulent conditions, one gets a saturation-like effect which 
is similar to that noted by Gracheva and Gurvich,® 4 and Dietz. This is 


‘See Appendix C. 
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Figure 18. Percent deviation of log-amplitude variance measurements 
vs skewness of log amplitude. 

particularly interesting since Eq. 95 was used to evaluate log-amplitude variance 
in References 91 and 95. Fig. 19 presents the difference a E 2 - cr 2 as a function 
of the third central moment. As pointed out by Hufnagel," ct £ 2 - a 2 = 2m 3 + higher - 
order terms. 

This relationship seems to be in fairly good agreement with the data. It is 
noteworthy that cumulative probability plots of log amplitude for various S values 
do not appear significantly different on Gaussian probability paper. Fig. 20 pre- 
sents results for typical data with S ranging from -0.143 to -0.694. 

Concluding this section, it appears that the validity of Eq. 95 is extremely 
sensitive to deviations from log normality and, in most interesting experiments, 
its accuracy is very dubious. 
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- m 3 x 10 3 

Figure 19. Difference cr £ 2 - cr* as a function 
of third central moment, m 3 . 

b. Saturation Level of Irradiance Fluctuations . In order to give more in- 
sight into the implications of Fig. 17, more data were recorded over the same 
range during strong turbulence conditions. An argon laser radiating at 4880 A 
was used. The output beam had a diameter of 1.2 mm between l/e 2 points and 
a divergence of 0.7 mr. The received energy was measured by a photomultiplier 
through a 1 . 5 -mm aperture. 

The results are shown in Fig. 21, where the circled points represent the 
argon data and the solid points are the He-Ne measurements already presented 
in Fig. 17. It is clear that cr T 2 fails to increase much beyond 0.3, regardless of 
how large cr E 2 becomes. At the time of this writing, there have been five 
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-0.143 


0.226 

-0.694. 


Figure 20. Cumulative probability of log amplitude 


experimental efforts known to the author directed toward measuring the satura- 
tion level for irradiance fluctuations at visible wavelengths. The net result seems 
to be that there is complete agreement that saturation occurs but very little 
agreement as to the actual level at which this happens. 

The work of Gracheva and Gurvich in the Soviet Union produced the first 
experimental evidence of the saturation of log-amplitude variance. The saturation 
level was measured as approximately 0.58. However, as pointed out by deWolf, 65 
there was an error in the data reduction, and the actual measured value was 0.14. 
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Figure 21. Comparison of log-amplitude variance :<r^ calculated from irradiance 
moments [Eq. 95] and a E 2 calculated from log-amplitude moments. The 

circled data points represent argon data, and the solid points are He-Ne meas- 
urements. The five levels shown are the saturation levels measured in five dif- 
ferent experiments. 


Following this, additional data were published by Gracheva and the saturation 

91 

level was observed to be 0.25. In the United States, Ochs and Lawrence pub- 
lished data indicating saturation at about 0.6. Two data -processing techniques 
were used; the first measured the slope of the cumulative density of log ampli- 
tude ; the second was entirely analog and used logarithmic amplifiers and RMS 
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voltmeters; the techniques gave similar results. Unfortunately, the authors 
concluded that their data were in agreement with those of Gracheva and Gurvich. 
Apparently they had not noted deWolf's correction to these data, although it had 

g 1 

been published nine months earlier. Later, using photographic techniques , Dietz 
measured saturation at about 0.21. In this case, as well as in Ref. 84, the var- 
iance of log-amplitude was computed by measuring irradiance statistics and 
transforming to log amplitude statistics using Eq. 95. The remaining experimental 
effort was conducted by Fried^ who observed saturation for both visible and in- 
frared fluctuations, and in addition noted that the saturation level was wavelength - 
dependent. The analysis technique based on the slope of the cumulative probabil- 
ity function was used here and resulted in a saturation level in the visible of 0.55. 
Since the experimental conditions in each of the listed instances were approxi- 
mately the same, it seems likely that only Fried and Lawrence and Ochs have 
really measured saturation. The saturation observed in the other experiments 
appears to be either a "statistical saturation" (i.e.,the use of Eq. 95) or perhaps 
an effect due to dynamic range limitations.* Dietz®® remarks that when his data 
were analyzed with the cumulative probability technique (which is less dependent 
on dynamic range), the values for log-amplitude variance became "somewhat 


*In addition to the published experimental efforts, there is a widely-quoted summarization of con- 
temporary propagation work by Strohbehn^ which has also contributed to the confusion. The 
saturation level is discussed and the factor of two error in the Gracheva/ Gurvich data is noted. 
The data of Gracheva are presented; through a succession of notation errors, missing factors of 
two and interchanging of variances (X^) with standard deviations (cTj ^ Strohbehn 

concludes that the saturation level is approximately 0.8 (p. 1307). It appears to this writer that 
the calculation should be = (0.8/2) = 0. 16. 
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higher," but not enough information is given so that meaningful comparisons can 
be made. 

In addition, there are some fundamental problems associated with the use 
of pulsed sources and photographic techniques such as those used in Ref. 89. Often 
the output of such sources is spatially quite non-homogeneous, and, as a result, 
the measured spatial probability density function of log amplitude is the convolu- 
tion of two densities: the output density associated with the source and the density 
associated with the atmospheric- induced fluctuations. In fact, the source charac- 
teristics may mask the atmospheric -caused fluctuation for ranges up to 200 
meters. 89 Other obvious problem areas include compensating for the nonlinear 
transfer functions of the film and overcoming the severe dynamic-range problem. 

Judged by the noticeable nonlinearity of the cumulative probability plots of log 

92-94 

amplitude presented in Ref. 89, as well as in other reports by the same 

group of experimenters, it does not appear that the combination of pulsed sources 
and photographic detection techniques permits accurate measurement of the 
density function of irradiance fluctuations. 

The difference o E 2 - a T 2 for the high -turbulence argon data presented in Fig. 

21 is superimposed on the original He-Ne data and presented in Fig. 22 as a 
function of m 3 . The data appear to be in approximate agreement with theory. 

We conclude that much of the data on the saturation of the variance of 
log amplitude are suspect.* The level at which saturation occurs appears 

2 

*The absolute accuracy of the very high-turbulence argon data (i.e. <7g > 0.3) presented in Figs. 

21 and 22 must be estimated with caution, because of the extended instrumental dynamic range 
required for such measurements (see Appendix C). However, the conclusions presented are not 
affected even if all data points with rtg >0.3 are deleted. 


81 







2.4 Spatial and Spectral Correlation of Irradiance Fluctuations 


This section presents the results of an experimental investigation of certain 
second-order statistics of the amplitude field. Specifically, we consider the 
two-point statistic (spatial correlation) and the two -wavelength statistic (spectral 
correlation) . Part (a) deals with spatial correlation, and Part (b) deals with spectral 
correlation. 

a. Spatial Correlation: Theory . We consider an electromagnetic wave 
which propagates through a medium characterized by zero conductivity, con- 
stant permeability, and a permittivity which is a random-point function. A 
degree of randomness is imposed on the amplitude and phase of the wave by the 
stochastic nature of the medium. The log -amplitude covariance for the wave as 
measured at points X and X 2 is 

(X,, x 2 ) = ^(Xj)}!^)- *6(x a )>}). (96) 

where 'C(X) is the log -amplitude value measured at X , and the angle brackets de- 
note an ensemble average over all the possible realizations of the permittivity 
field. When the electromagnetic field quantities are assumed to be locally 
homogeneous and isotropic , Eq. 96 becomes 

C ^(P) = ({^i) - (^X+P))}), (97) 

where 

P = I Xj - X 2 1 . 
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Tatarski 95 has evaluated Qg (p) for infinite plane -wave propagation through 
homogeneous isotropic turbulence, and Fried 9 ® has solved the spherical -wave 
case for similar turbulence conditions. Unfortunately, direct experimental test- 
ing of these results is difficult; existing optical detectors are square-law de- 
vices and, therefore, have outputs which are proportional to the square of the elec- 
tric field rather than the logarithm of the field. High-speed logarithmic processors 
with sufficient accuracy, moreover, are difficult to obtain and have, up to this time, 
limited the accuracy with which theoretical-experimental comparisons can be 
made. It is possible, however, to transform the results for log-amplitude covar- 
iance into irradiance covariance. This permits more direct comparisons of 
experiment with theory. 

The starting point is the covariance of irradiance Cj O) , defined as 

C^) = \{i(x 1 )-I 0 (x 1 )}{i(X 2 )-I 0 (X 2 )}^, (98) 

where 

M^i) = (ifxj). 

I 0 (x 2 ) = <l(x 2 )>, 

P = I Xj - x 2 1 , 

and the assumptions of homogeneity and isotropy have been incorporated in Eq. 98. 
From the definition of log amplitude -t(X) , we obtain 

I(X) =I 0 (X)exp[2^(X)] . (99) 
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We assume statistical homogeneity in the plane of the detectors ^I 0 (x i ) - I 0 ( x 2 ) 
s I 0 ] and Gaussian statistics for log amplitude. Then, 


C r iP) 


I 0 2 ^{exp [&,(*)] - l}{exp [ 2 1{\ 2 )] - 1 ^ 

I o 2 ^ ex P K,) + 24(X 2 )] - exp H* 2 )] - exp [2^(xj] + l^ • (100) 


The moment-generating function with parameters t x , 
distribution with random variables a 1 , a 2 is 

m ( t l' fc 2 ) “ ( eX P [tl +t 2 a 2 ]) 

= expftj + t 2 M 2 + \ (tj 2 a x 2 + 2rtj 


1 2 for the bivariate normal 


t 2 CT l 



(101) 


where 


( a i ) > 


<( a i“^i) 2 > 


( a 2 ) P'2 ’ 


<( C 


M 


! ) 2 > 


and 


<( a i ~P'l){ a 2 -P 2 )) 
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Therefore, the first term of Eq. 100 is the bivariate normal moment-generating 
function (m.g.f.) evaluated at t x = t 2 = 2, while the second and third are the 
m.g.f. for the corresponding marginal density functions. Using Eq. 99 and noting 
that 86 

<I(X)> = I 0 (X) = <I 0 (X) exp [24(X)]> (102) 


implies that 


<exp[2^(X)]> = 1 


and 

<t(X)> = - <^[^(X)- <^(X)>] 2 ^> - (103) 

Then, we obtain 

C t (P) = I Q 2 {exp 2/x 2 + + 2C^(p) +cr 2 2 )] - l\ 

3 1 o { ex P [ 4c -tO)] _ 1 } ’ (1 ° 4 ^ 

where C.£ (p) is the covariance function for log amplitude, and the p and a 2 symbols 
are respectively log-amplitude means and variances at X t and X 2 . The irradiance 
correlation coefficient y at separation p is 


Cj (p) _ exp [4C^ (p)] - 1 

exp [4C^ (0)] - 1 


(105) 
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and we define the log-amplitude correlation coefficient F(p) through 


Qe(p)/Qc< o) = f(p) . 


(106) 


Therefore, 


7(P) 


exp [4C^(0) F(p)] ~ 1 
exp [4C^(0)]-l 


(107) 


It is clear that 7(p) depends upon F(p) and also on the strength of turbulence 
through the scintillation statistic (0). Expanding the numerator and denom- 
inator of Eq. 107 in a Maclaurin series reveals clearly that in the limiting 
case, as the turbulence becomes very weak, the normalized log-amplitude and 
irradiance correlation coefficients become equal. This is intuitively satisfying, 
since for weakly fluctuating phenomena the logarithmic operator approaches 
linearity. 

Tatarski®® presents F(p) (Fig. 7-13) for plane-wave propagation through 
homogeneous and isotropic turbulence, and Fried derives a similar result for 
the spherical -wave case in Fig. 1, Ref. 96. We can now transform from log- 
amplitude correlation coefficient to irradiance correlation coefficient through 
Eq. 107 for any specified degree of turbulence [i.e. any value of (0)] . This 
has been carried out in Fig. 23 and Fig. 24 for the plane-wave and the spherical- 
wave case, respectively, and with log-amplitude variance values of 0.076 and 
0.460. These particular values were chosen to correspond with the experimental 
data, and, in addition, the 0.460 value is on the order of the upper limit of Rytov 
theory. The detector separation has been normalized by (kz) 1/2 , where \ is the 
optical wavelength and z is the propagation path length. 
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Figure 23. Irradiance correlation coefficient y as a function of normalized de- 
tector separation p/ (\z) 1/2 . The curves are obtained from Tatarski's 2 results 
for the plane -wave log-amplitude correlation coefficient. 


Experiment . The system for measuring spatial correlation is shown sche- 
matically in Fig. 25. A vertically polarized argon laser beam radiating at 

o 

4880 A with a divergence of 0.7 mr and a diameter of 1.2 mm was transmitted 
over a folded 1.17-km horizontal path. After reflection from a 51-cm flat mirror 
in the field, the beam was detected back at the transmitter station by two photo- 
multiplier tubes which had 1-mm entrance pupils and which were spaced from 
2 to 52 mm apart. The outputs of the detectors were amplified and then recorded 
on magnetic tape through an FM system with 2. 5-kHz bandwidth. For each de- 
tector separation, approximately two minutes of data were taken. 

Analysis of two-point statistics is complicated by the fact that each detector 
(and subsequent electronic chain) may have different scale factors. However, if 
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Figure 24. Irradiance correlation coefficient y as a function of normalized de- 
tector separation p/(Az) 1/2 . The curves are obtained from Fried's 96 results 
for the spherical-wave log-amplitude correlation coefficient. 


all gains remain constant during the time interval corresponding to a single data 
point (2 minutes), then gain variation can be easily taken into account. 

When the recorded data are played back, the voltages Vj and V 2 appearing at 
tape recorder outputs 1 and 2 are related to the irradiance values at detectors 1 
and 2 as follows: 


and 


V i = I x 

v 2 - k 2 i 2 , 


(108) 


where K x and K 2 include all detector and tape-recorder gains. Consider the 
playback of the data with outputs V x and V 2 connected to amplifier lines 1 and 2 , 
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respectively. Then the output voltage V 0 after a time interval equal to RC (the 
time constant of the averaging network) is 

v 0 = <A 1 A 2 A 3 (V 1 -V 1 )(V 2 - v 2 )> - V 01J (109) 

where and V 2 are the respective average values of Vj and V 2 , and A x , A 2 , A 3 
are amplifier gains. Similarly, if we repeat the playback with V t connected to 
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both amplifier lines, we get 


v 0 = <a 1 a 2 a 3 (v 1 -v 1 ) 2 > - v o: 


Repeating with V 2 connected to both amplifier lines, we get 


( 110 ) 


v 0 = <a,a 2 a 3 (v 2 -v 2 ) 2 > 


0 3* 


(HD 


Therefore , 


^ 0 l /(^02 ^03 ) 


1/2 = 


<(Vi- v 1 )(v 2 -v 2 )> 

{<(v 1 -v 1 )»><(v,-v l) »5r 


( 112 ) 


Applying Eq. 108, we get 


V 0l/( V 02 V 03) 1 '' 


2 = 


C, (P) 


{c Ii (0)c l 2 (0)} 1/ 


(113) 


Using the assumption of isotropy in the detector plane, we get 


c. (0) = C T (0) - C T (0). 


Therefore , 


V 


01 


{ v o2 v 03 ; 


\ 1/2 


y 


(114) 


So it is clear that one can evaluate y for each particular separation p by playing 
back the recorded data three times and measuring the appropriate voltages. 

From the two minutes of data available for each detector separation, a 30 -second 
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irradiance correlation coefficient 


portion was selected such that the signal levels were best contained within the 
linear range of the amplifier and multiplier circuits. Even so, it was apparent 
that some of the analyzed data exhibited peak irradiances which exceeded the 
dynamic range of the tape recorder. This problem will be discussed more fully 
in the next section. 

Experimental Results . Typical data are compared to spherical -wave theory 
and plotted in Fig. 26. Actual values for the log -amplitude variance at each data 
point were not available in this case, but the theoretical curves of 0.076 and 
0.460 can be expected to contain most of the data. For cases where the log- 
amplitude variance drops below 0.076, the pertinent theoretical curve is still 
very close to that shown for (0) = 0.076. If (0) becomes much greater than 
0.460, Rytov theory appears to lose its validity, and hence the analytical approaches 
which produced the results shown in Figs. 23 and 24 are of questionable usefulness. 



Figure 26. Spatial correlation of irradiance fluctuations. 
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The data in Fig. 26 show a reasonably good agreement with theory; how 
ever, there are some rather disturbing aspects. For most detector separations 
the data spread is rather large and is certainly greater than that which could be 
attributed to statistical sampling theory. Also, significant negative correlations 
for p/ (\z ) 1/2 > l appear consistently in the data of Fig. 26. 

Certainly some of the data scatter can be attributed to the lack of stationarity 
of the atmospheric statistics. It is well known that the refractive index structure 
constant can change by an order of magnitude within a few minutes of time. How- 
ever, this does not account for the data which lie above the 0.076 curve of Fig. 

26 because this curve is not significantly different from the curve for vanishingly 
small (0) obtained from Eq. 107. 

Another area of concern is the effect that atmospheric deviations from the 
Kolmogoroff model spectrum might produce on the measured correlation func- 
tion. It is shown in Appendix B that , if the refractive index structure function 
increases as r p (r being the distance between measurement points), then the 
log-amplitude variance takes on the range dependence L (3/2+p/2) and a wave- 
length dependence k (3/2_p/2 ) . Thus if p differs slightly from its usual value 
2/3,* then it is clear that the range and wavelength dependence of (0) are not 
drastically altered. However, it is not clear how seriously the shape of the co- 
variance function would change in response to varying p. 


’Lawrence and Ochs have suggested in private communications that near ground the refractive 
index structure function may frequently vary as r 1 ^ 2 due to the fact that input energy to the 
inertial subrange may be applied throughout the inertial subrange and not just at the outer scale 
as is usually assumed. 
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Additional data are shown in Fig. 27. In this case the data were processed 
through a digital computer (Raytheon 520) as well as through the analog instru- 
mentation of Fig. 25. The computer evaluated C,j> (0)by calculating the moments 
of the random variable 1/2 In [i/ ( I )], where I is the instantaneous irradiance 
at one of the detectors shown in Fig. 25 and ( I ) is the average irradiance at the 
same location. The value of log -amplitude variance changed from point to point 
for both the afternoon and the night data, but excursions from the mean value were 
nearly always less than an octave. Therefore, since the dependence of y on 
C£ (0) is rather weak, we can meaningfully form an "average log-amplitude 
variance" for both afternoon and night data. Some of the raw data were inspected 
in terms of their probability density function also. A typical example for one of 
the 30 -second data runs is shown in Fig. 28. The approximately linear behavior 
on log-gaussian probability paper indicates that the data satisfies this test for 
log-normality. 



Figure 27. Spatial correlation of irradiance fluctuations. 
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Figure 28. Cumulative probability of irradiance fluctuations. The data are plotted 

on log-normal probability paper. 


The results in Fig. 27 continue to exhibit the wide data spread noted previously. 
The nighttime data also show a significant negative correlation for p/(X z) 1 ' 2 
1.4. It is interesting that the strong turbulence data tend to take on a different 
shape which is characterized by increasingly positive correlation for large de- 
tector separations. This was observed several times in our data and has also 
been observed by others, either directly 9 ^’ 9 ® or indirectly. 99 In Ref. 97 it was 
found that y(p) evaluated at p - (\z) 1/2 may be as great as 0.40. 

It is tempting to conclude that this change in the measured shape of y(p) for 
strong irradiance fluctuations is real and is explainable on the basis of the break- 
down of Rytov theory. This may be the case, since it is generally agreed that 
Rytov theory does not successfully predict the behavior of (0) for very strong 
fluctuations; so one would not expect (0) y(p) to be accurate. However, 


95 


iiiiiiiiiiiiiiiniiniiiiiiiiiiiiii inn 


one must be very cautious about the results of covariance measurements made 
on such strongly fluctuating signals. Extremely large dynamic ranges are re- 
quired for the instrumentation,* and it is not always clear in published reports 
what the actual dynamic range capability is. 

From inspection of the experimental results shown in Figs. 26 and 27 and from 
the large data spread as well as the anomalous negative correlations which fre- 
quently occurred for p/('kz ) 1/2 > l, it was decided that there were three areas 
where the experimental technique could be improved. In an attempt to resolve 
whether the results had, in fact, been affected by the measurement technique, 
the following changes were implemented: 

(1) The dynamic range of the instrumentation was increased by 6 db. The 
limiting processor in the data reduction system was the FM magnetic tape re- 
corder. The signal representing instantaneous irradiance is inherently uni- 
polar, owing to the type of optical detectors employed. However, by DC -biasing 
the signal voltage to one extreme signal level, the useful recording range of the 
magnetic tape is fully utilized and is increased by 6 db over the unbiased case. 

The bias level must be subtracted during data reduction if accurate average 
values are required. The signal-to-noise ratio of the tape recorder in this con- 
figuration was 43 db. 

(2) The low-frequency cutoff implied by the capacitor coupling preceding 
the multiplier (see Fig. 25) was changed from 4 Hz to 1 Hz. 


*See Appendix C of this report and Appendix A of Ref. 46 for a discussion of dynamic range 
requirements in the measurement of log-normal random variables. 
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(3) The analog multiplier used in the data reduction (see Fig. 25) was re- 
placed by a digital multiplier. 

Steps 2 and 3 were actually implemented by replacing the preamplifiers 
(A lt A 2 ), the multiplier A 3 , and the RC averaging components with a digital 
processor: the Hewlett-Packard Model 3721 A correlator. The analog multiplier 
was a Burr- Brown Model 1671 and had a bandwidth of 40 kHz and a transfer func- 
tion accuracy of 0. 15$ . The bandwidth was not a constraint (the tape recorder 
was run at 7.5 inches per second and limited bandwidth from DC to 2.5 kHz) , but 
the transfer function accuracy was a source of concern. Especially for large 
p/(K z) 1/2 values, one can expect a significant portion of the covariance to be 
contributed when one input is at a very low level (where the percent error is 
high) and the other input is at a rather high level. Detailed specifications for 
the 3721 A correlator are not listed here, since they are readily available in the 
manufacturer's literature. However, we do note that it is a completely digital 
device with a low-frequency cutoff of either 1 Hz or DC and a high-frequency 
limitation of 250 kHz. It is possible to compute any of the following in an on-line 
configuration: 

(1) Autocorrelation of a single signal 

(2) Cross-correlation of two signals 

(3) Amplitude probability density function 

(4) Amplitude cumulative density function. 

All correlation calculations are made with 100 timelags of selectable spacing. 
Results are presented on a storage oscilloscope and can be interfaced with an 
X-Y recorder. 
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Additional data were recorded with a biased recording technique and processed 
through the digital correlator. The results are shown in Fig. 29. Log -amplitude - 
variance values were not available for these data, but both the morning and evening 
fluctuations were relatively weak and can be expected to lie closer to the 0.076 
curve than to the 0.460 curve. The quantity of data analyzed was small* but, even 
so, some significant differences are apparent. The data spread between runs 
(which were separated in time by 11 hours) are considerably less than those noted 



Figure 29. Spatial correlation of irradiance fluctuations (digital-data processing). 


*The digital processor was obtained on a short-term loan basis and was available only for two 
days. During this time the spatial correlation data presented in Fig. 7, as well as a somewhat 
larger quantity of spectral correlation data [to be presented in Part(b) of this section], were 
reduced and plotted. 
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previously in Figs. 26 and 27. Also, the anomalous negative correlations for 
P ^ (>vz) 1/2 are no longer present. Judged by this limited data sample, it ap- 
pears that the extension of the dynamic range and the employment of the seem- 
ingly more accurate digital processing have resulted in measurements which are 
in excellent agreement with spherical -wave theory. 

An added output which was available from this digital analysis was the more 
general time-dependent correlation coefficient y(p, r), defined as 

({!(*! - t+T )“ <I(X X , t +r))}{l(X 2> - t)- <l(X a , 

t + T )~ t + T )>} 2 ) ({i(x 2 . t) -<i(x 2 , 0>} 2 )] 

(115) 

where 

P = IXj-xJ • 

The assumptions of isotropy and homogeneity are again incorporated in the defini- 
tion; y(p, 0) is identical with the correlation coefficient y(p) defined in Eq. 107 
and used in Figs. 23, 24, 26, 27, and 29. 

Consider the following model for the propagation problem. The turbulons 
which exist along the optical path vary very slowly in time but are transported 
across the path in a horizontal sense by a wind vector V H N which is both hori- 
zontal and perpendicular to the optical path. * In this situation the field at X t at 


*We note that the component of the wind vector which is along the path has some effect on the 
problem, but it is shown in Ref. 2, p. 215, that this is usually insignificant. 
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time t j is the same as that at X 2 at time 


t 



where Xj and X 2 are position vectors which are collinear with one another and 
parallel to V H and which are also assumed to lie in a plane perpendicular to 
the optical path. The plus or minus ambiguity is easily resolved by inspecting 
the relative directions of X 1 - X 2 and V H N . If this model is accurate, one can 
always expect to obtain nearly perfect correlation between two side-by-side de- 
tectors (such as in Fig. 25) for some particular time delay 

lx x -x 2 | 


The time-dependent correlation function y(p, t) was evaluated for most of the 
data shown in Fig. 29. Some typical results are presented in Figs. 30 and 31. 

The data of Figs. 30 and 31. A were recorded about one hour after sunset. The 
wind was rather calm, with an average component of about 2.0 meters per second 
perpendicular to the optical path. This wind measurement was made at the 
transmitter-receiver station, and it is not known how well this measurement 
represents the actual average wind velocity along the path. It would have been 
more meaningful to measure the velocity at the folding mirror location, since 
turbulence in this region of the path can be expected to dominate the field 
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(B) 


Figure 30. Time -dependent spatial correlation 
coefficient y(p, r) as a function of time delay r. 

fluctuations, but this point was not accessible.* The Fig. 31. B data were recorded 
about one hour before sunrise, during a period of gusty winds with an average 
perpendicular component of about 7.5 meters per second. 


♦For spherical-wave propagation the turbulence near the source has little to do with the amplitude 
fluctuations, since in this region the beam cross section is small compared to the turbulence 
scale size. The net effect of this turbulence is to "steer” the beam as a whole rather than break 
up the beam. Near the receiver, the small-angle scattering aspect of optical propagation tends 
to minimize the effects of turbulence. 
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(A) 



(B) 

Figure 31. Time -dependent spatial correlation coefficient y(p, r) as a 

function of time delay r. 


The failure of any of the curves of Figs. 30 and 31 (or any of the other data) 
to exhibit near-unity correlation at some time delay shows clearly that the often- 
assumed propagation model described earlier is inadequate for predicting the 
time behavior of irradiance fluctuations. Several possibilities exist for explaining 
the observed data: 
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(1) The motion of the turbulons across the optical path may contain a signifi- 
cant vertical component due to buoyancy effects induced by the lapse rate. This 
effect can be expected to be emphasized in rather calm conditions (such as 
existed in Figs. 30 and 31. B). 

(2) The average wind component perpendicular to the propagation path may 
differ significantly along the path. 

(3) The concept of "frozen- in-turbulence" described earlier may not be 
valid. 

It is felt that (2) is probably not very significant, since for spherical-wave 
propagation it is the turbulence near the midpoint of the path which dominates 
and the average wind velocity need only be uniform in this region. Also, alter- 
native (3) seems unlikely since this concept (usually referred to as Taylor's 
hypothesis) has been extensively investigated by micrometeorologists 2 ^’ 25 and 
is generally considered to be a valid concept. We conclude that it is most likely 
the addition of a vertical component to the motion of the turbulons which has 
altered y(p, t) from its expected form. It is quite possible that greater degree 
of correlation would have appeared for some timelag if the detectors had been 
located such that the line joining their respective apertures was tilted rather 
than completely horizontal. We summarize the results of this section as follows: 

(1) It has been shown that the previously derived log-amplitude covariance 
can be transformed to the more useful irradiance covariance under the assump- 
tion of log-normality. 

(2) Measurements of the irradiance correlation coefficient over a 1. 17 -km 
horizontal path have been shown to be in good agreement with the Rytov spherical- 
wave theory. 
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(3) The propagation model employing "frozen-in-turbulence" undergoing 
horizontal transport has been found inadequate in this experiment under low- 
wind conditions. 

b. Spectral Correlation: Theory . We next consider the problem of evaluating 
the degree of correlation between amplitude fluctuations which are imposed on 
two beams of different frequency which have traversed exactly the same optical 
path. We neglect the dependence of the average index of refraction on wavelength, 
since the average value does not enter into the problem of amplitude fluctuations. 
Heuristically, different frequency beams will exhibit less than perfect correlation 
in their amplitude fields since the ratio of wavelength to turbulon size is different 
for each beam. Thus, even though each propagating beam "sees" the exact same 
refractive index field, the amplitude fields will tend to lose correlation if the 
measurement plane is far removed from a substantial amount of the turbulence. 
The spectral covariance of irradiance C t (0, 6) can be defined as 

c x (0 ,6) = (!yl 1 (X) - <I X (X)>| {l 2 (X)“ <I 2 (X)>}^>< (116) 

where 

I x (X) = the irradiance of beam 1 (of wave number k x ) , 

1 2 (X) = the irradiance of beam 2 (of wave number k 2 ) , 

and 


and we emphasize that both irradiances are measured at the same point X. As 
was noted earlier, for spherical- and plane-wave propagation, 

I(X) = I 0 (X) exp [24(X)] . (99) 
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(We drop the vector notation on the position vector X by applying the usual homo- 
geneity - isotropy assumptions.) Substituting Eq. 99 into Eq. 116, we write 

C r (0 , 9 ) = I Q1 I Q2 (exp [2^j + 24 2 ] - exp [2^J - exp [2^6 2 ] + l) • (117) 

Using the relationship 


(exp [2 4]) - 1 

and the right-hand side of Eq. 101, we have 

C T (0, 0) = I 01 I 02 [exp{4C^(0, 0)}- l], (118) 

where Cq> ( 0 , 6) is the log -amplitude covariance for beams 1 and 2 and is defined 
as 


C4(0, 0) = ({4 t (X)- (4 t (X)>}{4 2 (X)- (4 2 (X)>}> • (U9) 


In Ref. 100 it is shown that, for plane-wave propagation with wave numbers k x and 
k 2 through homogeneous and isotropic Kolmogorov-type turbulence, the generalized 
spectral covariance of log amplitude {p, 6) may be written as 


CliP, 0 ) = Ci ( 0 ) ( 1 - 0 2 ) 7/6 




5/6 , 


4 2 

1/2 

] r 4z e 

_k(l -0 2 )_ 

J 

lLk(i -e 2 )\ 


1/2 I 


( 120 ) 
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where F is the function use earlier in Eq. 106* and C i' (0) is evaluated 
at wave number |k 2 - k t |/k 2 + kj . Since F(0) = 1, (0, 6 ) becomes 

CiiO, 6) = Ci' (0) [l- e 2 ] 7/6 [l- 6> 5 / 6 ] 

" W(O)G(0) ■ (121) 

Substituting this result into Eq. 118, we have 

Cj(O,0) = I 01 I 02 [exp{4C^' (O)G(0)}-l] . (122) 

Evaluating Eq. 104 at p = 0, we have 

C x (0) = I 0 2 [exp{4C£(0)}- l] , 

and, substituting this into the defining equation for the spectral irradiance cor- 
relation coefficient y 12 (9) , we can write 

- c, ( 0 , 0)/{c lx ( 0 ) Cj 2 ( 0 )y /2 

exp (4C^' (0) G(0)} " 1 

^[exp [4C^(0)j - l] [exp {4C^(0)j - l]^ 2 (123) 

The result is now in a form such that for various strengths of turbulence (i.e. , 
various log-amplitude variances) the correlation between two plane waves with 

*We note that the function F(p) defined in Eq. 106 can be expressed more precisely in terms of 
the dimensionless variable p/(Az)^ 2 ; therefore, F (p) becomes F p/( Az) ^ 
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wavenumbers k x and k 2 can be numerically evaluated. We assume 7/6-power 
wave number dependence for log -amplitude variance (Eq. 52). Results are shown 

o 

in Fig. 32, where beam 1 is held constant at 6328 A and beam 2 varies over 
shorter wavelengths. The assumed turbulence levels apply to beam 1. The 
(0) values were chosen so as to bracket the conditions under which Rytov 
theory can be expected to be valid. The 0.05 (0) curve is not significantly 

different from the vanishingly weak turbulence condition. So now, just as in 
Part (a) , we have theoretical predictions which can be rigorously compared to 
experiment. 

Experiment . The experimental schematic is shown in Fig. 33. A helium- 

o 

neon laser operating at 6328 A was aligned collinear with an argon laser through 
a mirror-beamsplitter combination and transmitted over the folded 1.17 -km 
horizontal path. The two lasers had the following characteristics: 

Helium - Neon Laser (Spectra- Physics 130B) 


Power (TEM 00 ): 

Beam diameter (l/e 2 points): 
Beam divergence (total angle): 
Amplitude stability (short term): 
Amplitude stability (long term): 
Wavelength: 


1.0 mW. 

1.4 mm. 

0.7 mr . 

1% (RMS). 
2.8% (RMS). 
6328 A . 


Argon Laser (RCA - LD 2100) 
Power (TEM 00 , all lines): 125 mW. 

Beam diameter (l/e 2 points): 1.2 mm. 

Beam divergence (total angle): 0.7 mr . 
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Figure 32. Irradiance spectral correlation coefficient y l 2 (6). The assumed 
turbulence levels apply to beam 1, which has wavelength 6328 A in all cases. 
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Amplitude stability (total): 1.5% (RMS) . 

o 

Wavelengths: 4579 A individual power 5 mW. 

4658 A individual power 1 mW. 

o 

4727 A individual power 5 mW. 

4765 A individual power 16 mW. 

o 

4880 A individual power 50 mW. 

o 

4965 A individual power 10 mW. 

0 

5017 A individual power 5 mW. 

o 

5145 A individual power 30 mW. 

The argon unit was operated on one spectral line at a time. Each photomultiplier 
detector was equipped with optical filters such that detector 1 "saw" only the 
He-Ne wavelength and detector 2 "saw" only the particular argon wavelength 
being transmitted. The detected signals were amplified, as required, and recorded 
on magnetic tape in a manner identical to that described earlier. The two signals 
were later played back through the analog circuitry of Fig. 33 to evaluate 7 12 ( 0 )- 
Results . Data were recorded under conditions of strong (mid-afternoon) 
and weak (night) turbulence . The results are shown in Fig. 34 superimposed 
on the theoretical predictions of Fig. 32. Each plotted point represents 40 seconds 
of data.* Additional night data are shown in Fig. 35. The nighttime data show 
an acceptably small spread but exhibit very poor agreement with theory. The 


*We note that this part of the report actually compares plane -wave theory with spherical-wave 
experiments. As far as is known, no spherical-wave analysis of this problem is available. 
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Figure 33. Experimental schematic, spectral correlation. 

daytime measurements appear more compatible with theory but have a large 
spread. As before, no significant portion of this spread can be ascribed to sta- 
tistical sampling theory, since the 2B eq T product is large.* Both sets of meas- 
urements (Figs. 34 and 35) were obtained under weak -wind conditions. Typically, 
the wind speed ranged from absolute calm up to 2 meters per second. As has 
been noted by other experimenters, under such conditions there may not be 
enough turbulent energy present to generate an inertial subrange. Under such 

*See Appendix A for a discussion of the statistical aspects of the problem. 
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conditions, the turbulent field can be expected to lose its homogeneous and iso- 
tropic properties and take on something other than a Kolmogorov-type spatial 
spectrum. Solutions to the optical propagation problems do not exist for such 
conditions; therefore, measurements taken under such circumstances are 
nearly meaningless. Sufficient meteorological instrumentation was not avail- 
able at the time of measurement to determine definitely that this was the case ; 
therefore, this explanation is merely offered conjecturally. The conditions 
for the existence of an inertial subrange during midday are less dependent on 
wind speed, since the strong lapse rate which usually exists at this time is a 
potentially large source of input turbulent energy. It is possible that the after- 
noon data of Fig. 34 are more representative of a properly conducted experiment. 

After this preliminary look at data, it was decided that the remainder of 
the experimental measurements should not be taken at night; this was to minimize 
the possibility of encountering non-isotropic turbulence. In addition, the experi- 
mental modifications described in part a. were installed. That is, the tape- 
recorder dynamic range was increased by 6 db, the low-frequency cutoff in the 
covariance calculation was changed from 4 Hz to 1 Hz, and a digital processor 
was substituted for the analog multiplier and averaging circuits . Some range 
flexibility was also incorporated by installing another folding mirror on the range. 
This mirror was also flat, with a diameter of 30. 5 cm, and located so as to pro- 
vide a second folded path of length 400 meters. The results of measurements 
taken over this short path in late January are shown in Figs. 36-39. The graphs 
are labeled with the hours of the day of the respective data run, with zero hour 
at midnight. Each plotted point represents 30 to 60 seconds of raw data. A 
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Figure 37. Experimental results, spectral correlation— run # 2, digital 
data reduction. 
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Figure 39. Experimental results, spectral correlation— run # 4, digital 
data reduction. 
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typical set of measurements (i.e., a data run at each of the 8 argon laser wave- 
lengths) took about 30 minutes to complete. After each data point, the optical 
alignment between transmitted and received beams was checked. Both photo- 
multiplier detectors were checked at each wavelength to make certain that no 
cross-talk existed. 

The values for y 12 exhibit a very large scatter, although the large nighttime 

o 

correlations of Figs. 34 and 35 are absent. The 4727 A data point of Fig. 36 and 

o 

the 4965 A measurement of Figs. 37 and 41 were lost during the data analysis 
procedure owing to operator mishaps. 

The measurements taken over the full 1. 17 -km path are shown in Figs. 40-42. 
As in the short-range measurements the data scatter is so large that it masks 
any decrease of correlation with wavelength, if indeed there really is any over 

o o 

the span inspected (i.e., 5145 A to 4579 A). It is interesting to note that an order 
of magnitude change in the refractive index structure constant C n 2 would the- 
oretically be sufficient to move a data value from the C^(0) = 0.05 curve to a 
point somewhat below the (0) = 0.45 curve. In this light, the results obtained 
in Figs. 34-42 should not be too surprising. 

The average values of the short- and long-range data are plotted in Figs. 43 
and 44, respectively. The time -dependent correlation coefficient y 12 (6) defined 
as 


> 12 ( T - e ) 



> ( 124 ) 
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Figure 41. Experimental results, spectral correlation— run # 6 , digital 
data reduction. 
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Figure 42. Experimental results, spectral correlation— run # 7, digital 
data reduction. 
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Figure 44. Mean experimental results, 1. 17-km range. 
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was also evaluated for each data point in Figs. 36-42. In all cases y 12 (&, r) was 
maximum for r = 0. We summarize spectral correlation as follows: 

(1) The spectral correlation coefficient of log amplitude has been transformed 
to the irradiance domain in order to facilitate experiment-theory comparisons. 

(2) Measurement of spectral correlation coefficients over 0.4- and 1.17-km 
horizontal paths has not shown any decrease in correlation in the wavelength 
range 5145 A to 4579 A. 

o 

(3) A typical value for the spectral correlation coefficient between 6328 A 

o 

a nd 5000 A over a 1. 17-km horizontal path is about 0.7. 

2.5 Mean-Square Error and Chi-Square Testing of Irradiance Probability Density 

Functions 

The form of the density function which describes the amplitude fluctuations 
of an optical wave has recently been questioned on both theoretical 05 ' 101 and ex- 
perimental 101,102 grounds. Previously, owing to work in the Soviet Union 2 and 
the United States, 70,103 the log-normal distribution had been accepted as ac- 
curately describing amplitude statistics. Recent discussions in the litera- 
ture, however, have reflected a dissatisfaction with this result and have indicated 
that the Rayleigh or Rice distribution may be more accurate. In this section, the 
log-normal, Rayleigh, and Rice distributions are discussed and are compared with 
experimental data obtained with visible (0.488 fi) and infrared (10.6 /j) sources. 

The comparison between data and each of the three model distributions is car- 
ried out first through a mean-square error criterion (which emphasizes the 
modal region of the distribution) and secondly through a chi-square test (which 
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emphasizes the tails of the distribution). Each of the three distributions is dis- 
cussed below so as to establish its basis in optical propagation theory. 

a * Log-Normal Distribution . Following Tatarski (Ref . 2, Chap. 7), we con- 
sider a field E which is composed of an unscattered component 

E 0 = A 0 exp[jS 0 ] 

and scattered components 

E k = A k eX P[j S k]. 

where the index k ranges over a large number of independent contributors. 
Therefore , 



k 


and, for small fluctuations. 



Considering the real part of Eq. 125, we have 



E 


a : cos ( s k - s o) • 


(125) 


(126) 
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The left-hand side of Eq. 126 is the random variable log amplitude and, as the 
number of independent contributors (i.e., the range of k) becomes large, the 
Central Limit Theorem!*^ dictates that the amplitude statistics become log- 
normal.* 

Assuming Gaussian behavior for log amplitude, define its variance and mean 
as and a^, respectively, and note the following relationships: 


(a) 

. A 
In 

A o 

i 

_ T 

In r , 

A o 


(b) 

n = 

- 

> 

(127) 

(c) 

11 





Eqs. 127b and 127c follow from Eqs.103 and 104, with or ^ and I Q being the var- 
iance and mean of irradiance, respectively. 

Substituting Eqs. 127 into the Gaussian probability density function for log 
amplitude and transforming to irradiance quantities, we find the irradiance density 
g(I) to be 


g(I) 


i ( 87f a i ) 1/2 exp 


i 

8^7 



(128) 


*Note that the imaginary part of Eq. 125 is 


s - 


22 a; sin ( s u ~ s o)’ 


which implies a Gaussian behavior for the phase fluctuations when the range of k is large. 
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It is now clear that for log -normal -amplitude fluctuations the irradiance density 
can be evaluated by measuring the first two irradiance moments, solving for o 
through Eq. 127c, and substituting the result into Eq. 128. 

b. Rayleigh and R ice Distributions . -*' 05 In this case the total field is com- 
posed of a coherent component A 0 and an incoherent component which is repre- 
sented as a narrowband Gaussian process. That is, 

E(t) = [A 0 + x( t)] cosi) 0 t - y( t) sin a> 0 t , (129) 

where x and y are independent Gaussian processes with 

Elx] - Ely] = 0 

and variance cr 2 . Equation 129 can be rewritten as 

E(t) = r cos (o> 0 t + <£(t>) , (130) 


where 


2 4 - „2 


and 


<p - tan 


/( A o +x ) 


The probability density for r is called the Rice distribution and is written 
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Ri (r) 



(131) 
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where I 0 is the modified Bessel function of order zero. In the limiting case, as 
nearly all the energy is scattered, A 0 approaches zero and the probability density 
for r approaches the Rayleigh density Ra ( r) written as 


Ra(r) 



(132) 


Using deWolf’s^ results, when the field amplitude is distributed according to 
Eq. 131, we can write the irradiance density Ri' (r) as 


Ri' (r) 


1 




I - I ex p 
x 0 d 

L x o - i d J 

*0 

I 0 -*d 


(133) 


where 


I„ = E[I] 


and 



When the specular component becomes vanishingly small, I d - 0 and the irradiance 
density becomes 


Ra'(I) 


1 

"f exp 
x o 



(134) 


The irradiance densities specified by Eqs. 133 and 134 will be respectively called 
the modified Rice and the modified Rayleigh in the remainder of this discussion. 
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As before, it should be clear that if the amplitude field is assumed Bice- or 
Rayleigh-distributed the corresponding irradiance density can be easily gen- 
erated by measuring the first two moments of the irradiance fluctuations and 
substituting into Eq. 133 or Eq. 134. 

In summary, the approach used here is to measure the first two moments 
of the experimental data, and using these numbers generate the three model dis- 
tributions which would have been required in order to produce such experimental 
data. The three distributions are then tested to evaluate their compatibility with 
the measured density function. 

c. Probability Density Tests. Let the model log-normal, Rayleigh, and 
Rice distributions which are synthesized from the data be called g(I), Ra' (I), 
and Ri' (I), respectively. Then the mean-square error (MSE) of the fit is 

1 N 

n^H 11 )-^ 1 *)} 2 s mse ln , 

i=l 

N 

N ^{ d ( X i) - Ra ' (!i )} 2 = MSE Ra ', (135) 

i= 1 

and 

1 N 

n 2K d(ii ) ~ Ri ' (i) } 2 s “Hi" 

i - 1 

where d(l t ) is the experimentally measured density for the ith amplitude class 
and N is the total number of classes into which the fluctuating irradiance signal 
is partitioned. Clearly this test is not effective for evaluating the fit on the tails 
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of the distribution, since the quantity inside the bracket (Eq, 135) is always small 


in that region. 

In order to test the tails of the distributions, the chi-square goodness-of- 
fit parameter is employed. The defining equations are 

" { d ( I i ) ~ g l 1 i )} 2 
Xln * NAI 2_, g( I i) 


n {d(lJ-Ra' (ij } 2 

Xr V - nai ^ (136) 

and 



where Ai is the width of the amplitude class (assumed constant). 

The tails of the distribution are emphasized by normalizing inside the sum- 
mation by the model density. If the synthesized density is from the same family 
as the raw data, then it is well knownl04,106 that, as the sample size becomes 
large, x 2 takes on a chi-square density with N-l-S degrees of freedom, where S 
is the number of parameters of the model density which were estimated from the 
data sample. In this case, for any a the A can be evaluated which satisfies 

p k 2 -i-s >A ] = a • 

If we define the null hypothesis H 0 , that the target population had a density be- 
longing to the same family as the model distribution, then H 0 is rejected if Xn-i-s 
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is greater than A and is accepted if H 0 is less than A. The rejection of H 0 when 
it is actually true is called a Type I error and has probability a . The accepting 
of H 0 when it is actually false is called a Type II error, but nothing definite can 
be concluded about its probability unless a particular form is assigned to the 
target density. 

d. Experimental Analysis . Two types of experimental measurement will be 
analyzed in what follows. The only difference is in the type of optical source 
utilized. We first discuss the results obtained with a C0 2 laser radiating at 
10.6m and follow with the results obtained with an argon laser operating at 
0.488m . 

The 10.6m signal analysis is based on the same raw data that were used in 
Section 2.2 for log-amplitude variance analysis. However, it was apparent from 
the statistical analysis in Section 2.2 that longer time records would be required 
for density function estimation. A 40 -second record was the maximum practical 
length for computerized processing. The amount of data available was seriously 
reduced because of frequent shifts in the laser output power. All the data were 
inspected on a chart recorder, and only those 40 -second increments which were 
free from transmitter power fluctuations were analyzed. The detected irradiance 
signal (described in Section 2.2) was sampled at a 1-millisecond rate by an 
analog -to-digital converter and fed into a computer. The computer divided the 
amplitude range spanned by the signal into 200 evenly spaced classes. The prob- 
ability density function d(I) and the first two moments were computed. These 
moments were then substituted into Eqs. 128, 133, and 134 in order to generate 
g(I), Ri ' (I), and Ra' (I). The mean-square error calculation was then carried 
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by means of Eq. 135, and the chi-square calculation was performed by means of 
Eq, 136. The results are shown in Table XI. 


TABLE XI 

Mean-Square Error and Chi-Square Test Results of Weak- 
Scintillation Data (A. = 10.6 fj . ) 


Data 

Record 

Duration 

(sec) 

Mean-Square Error 

Chi-Square 

mse ln 

MSE Ra - 

MSE Ri - 

Xln 

^Ra' 

*Ri' 

#1 

40 

0.247 X 10 -1 

0.319 x 10 1 

0.211 x 10 -1 

0.243 x 10 3 

0.252 x 10 5 

0.139 x 10 3 

#2 

40 

0.175 X 10 _1 

0.260 x 10 1 

0.250 x 10 - 1 

0.192 x 10 3 

0.230 x 10 5 

0.228 x 10 3 

#3 

40 

0.899 x 10 -1 

0.111 x 10 2 

0.957 x 10' 1 

0.179 x 10 3 

0.252 x 10 5 

0.182 x 10 3 

#4 

40 

0.218 

0.168 x 10 2 

0.278 

0.231 x 10 3 

0.259 x 10 5 

0.291 x 10 3 

#5 

40 

0.184 

0 . 903 x 10 1 

0.127 

0.866 x 10 3 

0.255 x 10 5 

0.341 x 10 3 


In all cases the modified Rayleigh is a poor fit to the data. The log-normal and 
modified Rice appear nearly equally matched to the data from both a mean-square - 
error and a chi-square criterion. In Table XII, ratios of mean-square error and 
chi-square are shown for each record. Also shown is the measured scintillation 
statistic cr j/i o for the experimental data. The average value of the mean-square 
error ratio is 1.01 and the average value of the chi-square ratio is 1.38. 

The three model densities and the experimentally measured density for 
record #1 are shown in Fig. 45. Extra points are plotted near the mode where 
the data scatter is greatest. In places, the modified Rice is so close to the log- 
normal that they cannot be visually separated on a graph of this type. Table XIII 
is a portion of the computer printout for Record #2. Columns 6-8 are the 
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TABLE XII 


Mean-Square Error and Chi-Square Ratio for Weak -Scintillation Data 


Data Record 

(MSE ln /MSE r ,') 

UWi) 

a i /l o 

#1 

1.17 

1.75 

0.119 

#2 

0.70 

0.84 

0.126 

#3 

0.94 

0.98 

0.086 

#4 

0.78 

0.79 

0.086 

#5 

1.45 

2.54 

0.095 


cumulative distributions for the data, log-normal, and Rice densities. Columns 9 
and 10 list the contributions to x L n and x R 2 ' from each of the 200 classes. 

It is apparent that, under mild fluctuation conditions (cr /l 0 ~ 0 . 1) , the log- 
normal and Rice densities are very similar and are equally suitable for describing 
amplitude fluctuations. 

The visible data were obtained over the 1.17-km path with the argon laser 
(X- = 0.488 fj.) described earlier. The detector had a 1.0 -mm aperture and pro- 
duced a signal which was recorded on magnetic tape through an FM module. The 
signal level was biased so as to maximize the tape-recorder dynamic range (as 
described in Section 2.4a). The bias was subtracted from the data during the 
computer analysis. Essentially the same analysis was performed on these data 
as that just described for the C0 2 data. However, the interpretation of the re- 
sults is more complicated, owing to problems caused by dynamic range limitation 
in the instrumentation. The fundamental problem is this: the signal which is 
analyzed by the computer is actually the sum of a noise signal (tape recorder 
noise, photomultiplier noise, etc.) and the desired irradiance signal. Considering 
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Figure 45. Synthesized log-normal and modified Rayleigh and Rice distributions 
for weak irradiance fluctuations. The normalized RMS irradiance fading 
(o'j/Iq ) is 0.119. 
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Computer Results for Record #2 of the Weak-Scintillation Data 
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the noise and irradiance signals independent, we conclude that the probability 
density of the composite signal is a convolution of the individual densities as- 
sociated with the noise and irradiance signals . When the irradiance fluctuations 
are strong, there is a significant amount of time when the irradiance value is near 
zero. When the noise is summed with these near-zero values, the resultant fre- 
quently takes on negative values. In order to compensate for these physically 
impossible "negative irradiance" values, the d(l . ) , i = 1, 2, . . . 200, as printed 
out by the computer, was shifted so that d(l 1 ) was assigned to the first positive 
class, d(l 2 ) was assigned to the second positive class, etc. 

An additional problem area was that the tape recorder noise was not a con- 
stant but depended on the signal level being recorded. A calibration sequence 
consisting of a plus, zero, and minus DC voltage was recorded on the magnetic 
tape just prior to each data run. During data analysis, the variance associated 
with each of these levels was analyzed so as to evaluate the tape-recorder noise 
at each record level. In this experiment the negative calibration was very near 
the zero irradiance level for the system, and therefore it was this noise level 
which was important. For chi-square calculations, it was noted that the first 
few classes of the irradiance density (which were seriously contaminated by the 
noise) were dominating the chi-square value for each of the three model densi- 
ties. For this reason, the first two or three classes (depending on the measured 
noise) were omitted from the chi-square calculation in all cases. 

The results of 23 data runs of 20-seconds length are shown in Table XIV. 

The standard deviation divided by the mean irradiance is also listed for each 
record. Note that when o^/ig is greater than 1, it is not possible to generate 
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TABLE XIV 



#1 0 . 304 x 10' 2 0.249 0.178 

2 0.413 x 10 2 0.273 0.161 

3 0.152 x 10 1 0.401 0.295 

4 0. 116 x 10 1 0.358 0.269 

5 0 . 286 x 10“ 1 0.552 0.333 

6 0.263 x 10" 1 0.829 0.416 

7 0.240 x 10" 1 0.648 0.362 

8 0. 160 x 10" 1 0.527 0.276 

9 0.369 x 10" 1 0.439 0.420 

10 0.266 xlO" 1 0.389 


4,180 0.885 

10,810 0.892 

195,000 0.833 


0.800 


4,743 17,300 0.823 

4,480 625,000 0.812 

4,800 6,840 0.970 


1.005 


11 0.851 x 10" 2 0.554 0.257 

12 0.113 x 10" 1 0.697 0.238 

13 0.440 x 10" 2 0.447 0.158 

14 0.400 xlO" 2 0.426 0.191 

15 0.220 xlO" 1 0.378 0.310 

6 0.209 x 10"' 0.322 



0.806 


13,100 0.732 

4,118 0.762 


5,222 


0.915 


17 0.203 x 10"' 0.440 0.273 

18 0.740 x 10 2 0.378 0.227 

19 0.220 xlO" 1 0.360 0.319 

20 0. 105 x 10" 1 0.354 0.307 


9.U0U 0.838 
15,540 0.840 
4,810 0.940 


0.935 


21 0. 190 x io i 0.394 0.342 


0.932 


22 0.118 x 10 1 0.416 0.279 


0.865 


23 0.128 x io i 0.710 0.369 


2,100 0.815 
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a Rice density using deWolf's results (see Eq. 133). The results show that in 
all cases the log-normal density is a better fit to the data from both a mean- 
square-error and chi-square criterion. Fig. 46 plots the ratio of mean-square 
errors for the modified Rayleigh and log-normal distributions versus strength of 
the irradiance fluctuations (o-j/T 0 ) . The ratio decreases substantially as the 
fluctuations grow stronger, but there is still an order -of -magnitude difference 
for very strong fluctuations. Fig. 47 plots the MSE ratio for the modified Rice 
and log-normal distributions and shows the inadequacy c£ the former. The 
value of the log-normal MSE is shown in Fig. 48, and no clear trend is evident. 
Corresponding plots for the three chi-square values are shown in Figs. 49-51. 
These graphs also show clearly that the log-normal distribution is a much better 
fit to the data than the other models used. The three synthesized densities and 
the data of data record 5 (Table XIV) in Fig. 52. Every class value of the data is 
plotted up to I =0.2, and above this value only every fifth class is shown (for 
clarity purposes). A typical computer printout (data record #3) is shown in 
Table XV . The column headings are identical with those of Table XHI except 
for 8. and 10. , which now pertain to the modified Rayleigh density. The first 
two classes were omitted from the chi-square calculation, owing to the noise 
problem discussed previously. 

Neither the C0 2 data nor the argon data have been discussed up to this point 
from the standpoint of hypothesis testing. In Table XVI are listed the percentage 
points of the chi-square distribution for 200 degrees of freedom.* So it should be 


•The A values for 198 or 199 degrees of freedom are not significantly different from those listed 
in Table XVI. 
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NORMALIZED RMS FADING, ^ 


Figure 49. Chi-square ratio for the modified Rayleigh and log-normal 
distributions as a function of RMS fading,^/!,, . 
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Table XV. Computer Results for Record #3 of the Strong -Scintillation Data 

Mean-square error and chi-square test results for strong 

scintillation data (X = 0.488^). w 
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(Continued) 
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Table XV. (Continued) 
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TABLE XVI 


Percentage Points of the Chi-Square Distribution 
With 200 Degrees of Freedom 


^[*200 > A] a. 

a 

A 

0.20 

217 

0.10 

226 

0.05 

234 

0.025 

241 

0.01 

249 

0.005 

255 

0.001 

267 


clear that for every data record tabulated in this section, the hypothesis that the 
data sample came from a log-normal, Rayleigh, or Rice distribution can be tested 
at any level of significance (i.e. , any a value) . However, it is doubtful that this is 
a very useful approach in this experiment since in each data record a large num- 
ber of classes had zei’o observations. One of the generally accepted rules for 
applying hypothesis testing to the chi-square goodness -of -fit parameter is that 
every class of the measured density function contain one or more observations. 

In this experiment we attempted to comply with this requirement by coding the 
analysis so that the first class contained the lowest irradiance signal en- 
countered in the data and the last class contained the highest irradiance level 
in the data. Even so, there were sometimes as many as 80 and often about 50 
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classes which had zero observations at the end of a 20-second data run. So 
it is doubtful that hypothesis testing has much utility under these conditions. 

This section has shown that the probability density function for irradiance 
fluctuations has been tested under very weak and strong fluctuation conditions. 
For weak fluctuations, the log-normal and Rice distributions were found to be 
equally compatible with the data; for strong fluctuations, the log-normal dis- 
tribution appears to fit the data much more accurately than either the Rice or 
the Rayleigh distribution. 

2.6 Summary and Conclusions 

This dissertation has reviewed the optical propagation problem from first 
principles and has presented a listing of available theoretical results for various 
source types and propagation paths. The Obukohov-Kolmogorov model of tur- 
bulence has been emphasized, since it appears to be the most realistic and ac- 
curate statistical description. 

The results of a group of experiments conducted over 0.4-km and 1.17-km 
near-ground horizontal ranges have also been presented. The experiments were 
directed toward measuring parameters which control the quality of a communica- 
tion channel. Results were compared to Rytov-based propagation theory as de- 
veloped in the United States and the Soviet Union. 

Signal fading characteristics as denoted by the variance of log amplitude 
were measured simultaneously for visible and infrared beams which had tra- 
versed the same path. The results were in excellent agreement with theory and 
show clearly the reduction in fading which can be obtained with longer -wavelength 
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signal sources. The valid parameter ranges of Rytov theory are known to be 
limited by the saturation of the log-amplitude variance. Published results for 
the saturation level are contradictory. Early experiments indicated saturation 
at about 0.2, but recent results indicate that the saturation level is closer to 0.6. 
A possible explanation suggested by our experimental data for this discrepancy 
has been discussed. Data reduction techniques which assume log-normality for 
amplitude statistics appear to be very limited in accuracy. 

Diversity communication techniques have been investigated through the 
spatial and spectral irradiance correlation functions. The respective log-ampli- 
tude correlation functions were transformed to the irradiance domain so as to 
permit direct experiment-theory comparisons. The spatial correlation measure- 
ments showed good agreement with theory, but an extreme sensitivity to instru- 
mentation accuracy and dynamic range was noted. Judged by these measurements, 
it appears that the correlation distance for irradiance fluctuations can be ac- 
curately evaluated through Rytov theory. The spectral correlation measurements 
exhibited an unacceptably large scatter and did not show the wavelength dependence 
predicted by theory. A typical measured correlation for sources at 0.633m and 
0.488m which had traversed the same path was 0.7. This indicates that spectral 
diversity techniques in the visible will not be an effective means of overcoming 
turbulence-induced signal fading. We note, however, that the same measurements 
indicate that the coherent bandwidth limitation which exists at radio frequencies 
is not a significant problem at optical frequencies. Therefore, atmospheric 
turbulence will not impose a limitation on the information bandwidth of an optical 
carrier from a coherency standpoint. 
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Analysis of the quality of a communication system (particularly a digital 
system) requires knowledge of the density function of the amplitude statistics. 
Data have been obtained, analyzed, and compared to three suggested model 
densities (log-normal, Rayleigh, and Rice). A mean-square error and a chi- 
square test have been applied as the criteria of comparison. Under weak fluctua- 

i 

tions, the log-normal and Rice distributions were found to be equally well suited 
to the data; under strong fluctuations, the log-normal distribution was an order of 
magnitude closer to the data in both the mean-square-error and the chi-square 
tests. 
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Appendix A 


Determination of the Number of independent Samples for Processes With 
Arbitrary Spectral Densities 


In this section we consider the problem of finding the number of independent 
samples N that one obtains in a continuous data run of length T taken from a 
random process x( t ). 

First consider the case where x has a one-sided spectral density G x ( f ) which 
is flat -baseband , i . e . , 


K 

G x ( f ) = g for 0 < f < B 

and 

G x (f) = 0, for f > B ; 


(Al) 


then, the autocorrelation of x, R x (r)ia 


* (-) 


■ r 

j 0 


G (f)cos277f T df - K 


sin 2n B t 
2ttBt 


(A2) 


So, 


. R x (0) = e[x 2 ] = K (A3) 

and the first zero of R x (t) occurs at r - 1/2B. Therefore, samples taken at 
times 1/2B apart are uncorrelated and (for normal processes) independent. In 
a data rim of length T, it is now clear that one can obtain at most T/1/2B = 2BT 
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independent samples. Note also that 


/•CO >*co 

1 1 sin 2ttBt K _ K 

T J R x ( T ) dr “ t I K 27rBr dr - 2BT - N • (A4) 

Now consider the case where G x (f) is not flat -baseband but, rather, has an 
arbitrary shape. Then one can evaluate R x (r) through the usual Fourier 
Transform relation. Then 


_1 

T 



R x ( r ) dr 


R x (°) 

2B e q . T 


(A 5) 


where B eq is the equivalent rectangular bandwidth of the process. Therefore, 


B 


R x (0) 


/•CO 

2 R x ( t) dr 



1 


■y(r) dr 


(A6) 


where y(r) is the normalized correlation coefficient defined for stationary 
process as 


E[x(t) x(t + T)] 

“ VAR (x) ■ (A 7) 

So in this case the number of independent samples N obtained in a record of length 
T is 


N 



y(r) dr 


(A8) 
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In practice, two approximate techniques permit a quick evaluation of the above 
integral. Many physical processes have autocorrelation functions which can be 
closely approximated by the exponential fit e~ a|T| . In this case 

aT 

N = ~2 ' (A9) 


Alternately one may approximate y(r) by a linear fit, i.e. , 


y(r) - 1 - a| r 


0 — | T | 



= 0 , 


(A10) 


elsewhere. In this case, 


N = aT . 


(All) 
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Appendix B 


Range and Wavelength Dependence of Log-Amplitude Variance for 
Arbitrary Turbulent Spectra 

Consider a refractive index structure function of the form 

D n (x) = a 2 x p , (Bl) 

where a 2 is an arbitrary constant and p is also a constant in the range 0 < p < 2. 
Using Tatarski's^ results, the one -dimensional spatial spectrum corresponding 
to Eq. Bl is 


w(/< ) 


r ( 1 + P) 
2rr 


. ^P 
sin — 2 ~ 


a 2 c~<P +1 > *Up + 1> 


(B2) 


For isotropic random fields, the three-dimensional spectrum [<£(*)] is related 
to the one -dimensional spectrum [w(k)] through 




1 dw( k ) 
2jtk ck 


(B3) 


Therefore, 


<£(*) = A/<' (p+3) , 


(B4) 


where 


_ r ( 1 * P) 
( 2 n) 2 


77 P 


, 2 n _ (P + U 


(P+ 1) 
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For the case where the field quantities are homogeneous and isotropic in planes 
perpendicular to the direction of propagation, the covariance of log amplitude 
can be written as 


Q(r) 



(«r) F(k, 0) k dx , 


(B5) 


where F(«, 0) is the two-dimensional spatial spectrum in the measurement plane. 
For infinite plane -wave propagation through homogeneous turbulence, Tatarski^ 
has shown in his Eq. 7.50 that 


F(k, 0) = 77k 2 L I 


K 2 L 


sin 


4>(k) » 


(B6) 


where 

k is the wave number of the optical carrier, 

and 

L is the path length through the turbulence. 

Substitute Eqs. B4 and B6 into Eq. B5; the log-amplitude covariance becomes 


Qfc(r) - 277^ A k L J o (Kr) 1 


r 

J a 


x 2 L 


sin 


'( p+2 > ck , 


(B7) 


and the variance is 


C*(0) 


2t 7 2 A k 2 L I f 1 - 
■' o ' 


k k 2 l\ 

K 2 L Sln k ) K 


< p+2 ) d< . 


(B8) 
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Let 


and 


£ 


K 2 L 
k 




2L 

k ck » 


and substitute into Eq. B8 . Then 


Qt(0) 


2n 2 A k 2 L 



1 - 


^)(i)(£r 3> % 


= 77 2 Ak 3/2 -p /2 L 3/2+p/2 (l- ^i)^-( p+3 )/2 . (B9 ) 

-'o 


So, for infinite plane-wave propogation, the log-amplitude variance exhibits a 
k 3/2-p/2 wave-number dependence and a L 3/2+p/2 range dependence. For spheri- 
cal-wave propagation through homogeneous turbulence, the log-amplitude variance 
equation takes the form 
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where C(x) and S(x) are the Fresnel integrals 


and 


C(x) 


S(x) 



Let 


£ 


k 2 L 
2rrk ’ 


and substitute Eq. B4 into Eq. BIO. Then 


ty(0) 


k (3/2-p/2) L (3/2+p/2) 




+ 



^-<P + 3)/2 . 


So, again, the log-amplitude variance has a wave -number dependence k 3/2 ~ p/2 
a range dependence l 3/2+p/2 . 


(Bll) 


and 
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Appendix C 


Dynamic Range Considerations in the Measurement of Log- 
Amplitude Variance 


One of the major problems confronting the experimenter as he investigates 
amplitude fluctuations at optical frequencies is the dynamic range of the instru- 
mentation. The amplitude fluctuations are characterized by short periods of 
very high level (compared to the mean) and somewhat longer periods of low signal 
level (compared to the mean). Accurate recording and meaningful statistical 
analysis of this type of phenomena require extended instrumental dynamic range. 
The problem has recently been looked at analytically for the first time in Ref. 46 
In this analysis, it was assumed that the log -amplitude variance o-/ was evaluated 
through irradiance statistics by means of 


a 


2 



(Cl) 


where 

Io = E(D 

and 

CT i 2 = E ([l-Io] 2 )- 

It was also assumed that the available dynamic range was evenly divided about I 0 . 
That is, if I H and I L are the highest and lowest signals that can be handled by the 
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instrument, then the fluctuating signal is sealed so that 



A numerical integration of the right-hand side of Eq. Cl was carried out for various 
<y and dynamic range values. The results are presented in a graph of measured 
log-amplitude variance ( cq> 2 ) versus actual log-amplitude variance ( o -/ ) for 
various dynamic ranges. A typical result is that an 



18 


is required to measure a relatively weak 0 / = 0.2. This range (18 2 ) is approx- 
imately the state-of-the-art in FM magnetic tape recording. 

It appears that that these results may be unduly discouraging. If one is actually 
measuring erg 2 through Eq. 1, then E(I) and E(l 2 ) are the only statistics of inter- 
est and the available dynamic range would not be alloted equally about I 0 ; rather, 
it would be biased so as to accommodate peak signals. In addition, the numerical 
integration which produced the results did not clip signals which exceeded the al- 
lowed range but, rather, set these signals equal to zero. The amount of error 
introduced by this technique could be substantial. In this appendix, a different 
analysis is presented which may be more accurate in some cases. 

It is assumed that the log-amplitude variance is measured not through Eq. Cl 
but, rather, by evaluation of the probability density function for the variable log 
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amplitude {•i i ) , where 


and 



= 1 , 2 , 


n = sample size , 


(C2) 


and by computing moments e[-£. ] , EjJ'Cj ~'£ i ) 2 J in the usual manner. In this 
case it is efficient to allot the dynamic range as above, i.e., 


I 

T 


H 

0 



R . 


Our model for the system is 



l. 

1 out 

, for 

Ui 1 < R , 

in 

' in^ 

l. 

^ In R = 

B, for 'l. > R , 

1 out 

1 in 

and 

- tj In R 


1 out 

= - B, for <i < 

in 


(C3) 


The system is shown in terms of irradiance levels in Fig. Cl. We assume log 
amplitude to be normal with mean /j-% and variance cr£ . That is, l i is n(^, a£ ) . 
So, the measured mean-square value of •i is 


(•b r~B /•“ 

^^In.eas. = ^ n ’ a l ) M + J B 2 n(^, ) d-t + B 2 n(/^, ol) di. 

J-B •'-CO •' R 


(C4) 


Let 


y - l ~ fj.% . 
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Then, yisn(0, cr{f), and 


E[l 2 } 



y + n(0, a£) dy 



(C5) 


Consider the second and third integrals. If we assume B > , then we can 

change the limits with small error so that 



n(°, cr^) dy 




n (°, cr£) dy | 



(C6) 


where 


erf (x) 



The first integral in Eq. C5 can be expanded and evaluated by parts. 



nY n {°' 


a l ) dy 



+ 2i±i y + H-l ) n (0 , erg ) dy . 


(C7) 
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The first part is 



where 


P(x) 



e~ t2/2 dt . 


The second part is 


rB-u-l 

I y n(0, cr£) dy 



+ 


r B ~n 

Jb+h,£ 


y n (O , a£) dy . 


(C9) 


The first integrand is odd and therefore vanishes. Therefore , 


f B '*e r B ~^l _ / B \ 

I 2^y n(0, CTj?) dy = 2fi£ yn(o, o^ 2 ) dy = " B z(— J . (CIO) 


where n (o , crj? ) in the interval ( B + to B - has been approximated by its 
value at the midpoint B, and 


Z(x) 


7^ exp 




The third part is 


r*-n 

J-b-^ 


n( 0 , erg) dy 



2 



0 , a^) dy . 


Let 





Then, 


r B ~n 

y 2 n ( 0 , cr£ ) dy = 



where P(x) is the complete Gamma Function 


r (x) = 



dt 


and I(u, p) is the incomplete Gamma Function Ratio 


I(u, p) 


! r 

f(p + 1 ) t p e’ 1 dt . 


(Cll) 


(C12) 
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Simplifying and combining terms, we have 




meas. 


= I (7^ •*) - W B z(l) 



(C13) 


Realizing that the dynamic range constraint is much less severe for measuring 
EM , we assume EM meas = then 




*£ I 


B J 


/6 




1 

1 


b z 



-B- 

a l 



+ B 2 



(C14) 


The result is now in a form such that various strengths of turbulence ) and dy- 
namic ranges (B) can be put in and the measured strength of turbulence (cr£ } 

\ meas./ 

evaluated from tables. The calculation of the right-hand side of Eq. C14 requires 
that one recognize that - - o-^, which has been shown in Eq. 103 of Section 2.4. 

As an example, consider a data analysis system which has a signal-to-noise 
ratio of 43 db. As an analog system, it has a definite upper limit ( I H ) but its 
lower limit is not defined. However, the noise level of the system may be defined 
as the smallest signal level which can be accurately analyzed. Thus, 


43 


20 log j— 


J-H ! 0 

20 log t~ 

A 0 A L 


40 log R . 
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So, 


R = 11.9, 


and 


- I» 2 

- 2 fa Iq 


2-tnR = 1 . 24 


Eq. C14 can now be evaluated for various <x£ 2 values. The results of this calcula- 
tion are shown in Fig. C 2, where the points represent the actual calculations . 

Some caution must be used when evaluating Eq. C14 for small B/o^ values, 
since the approximations used in the derivation may lose accuracy. For B /cr^ > 2 , 
the first term in Eq. C14 dominates. Specifically, if 




B _ meas. 

2 , then, 


= .75 ; 


if 


B 

- 2.5, then, - 




meas . 

- = .95 


a l 


Defining the dynamic range (DR) of a device as its signal -to-noise ratio, 51 


DR 


20 log r r 


20 log R 


2 


40 log R , 


(C15) 


B = 2 In R , 


*With this definition, our dynamic range is numerically twice that used in Ref. 46 with units in dB. 
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Therefore , 


DR = 40 log e 2B = 80 B log e 


(CIS) 


= 34.7 B . 


Considering two levels of desired accuracy (5%, 25%) in the measurement sys 
tern, then 


and 


DR - 86. 7 cr£ (5% Accuracy) I 

(C16) 

DR = 69.4 cr£ (25% Accuracy) . II 


Fig. C3 plots I and II. 




o 



Figure C3 


NASA-Lauigley, 1971 23 
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